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ABSTRACT 

We trace the evolution of central galaxies in three ~ 10^^ Mq galaxy groups simulated at high 
resolution in cosmological hydrodynamical simulations. 

In all three cases, the evolution in the group potential leads, at z = 0, to central galaxies that 
are massive, gas-poor early- type systems supported by stellar velocity dispersion and which resemble 
either elliptical or SO galaxies. Their z 2 — 2.5 main progenitors are massive (M* ~ 3 — 10 x lO"*^^ 
M0), star forming (20 — 6OM0/yr) galaxies which host substantial reservoirs of cold gas (~ 5 x 10^ 
Mq) in extended gas disks. Our simulations thus show that star forming galaxies observed at z ~ 2 
are likely the main progenitors of central galaxies in galaxy groups at z = 0. 

At z ~ 2 the stellar component of all galaxies is compact, with a half-mass radius < 1 kpc. 
The central stellar density stays approximatively constant from such early epochs down to z = 0. 
Instead, the galaxies grow inside-out, by acquiring a stellar envelope outside the innermost ~ 2 kpc. 
Consequently the density within the effective radius decreases by up to two orders of magnitude. Both 
major and minor mergers contribute to most (70l^5%) of the mass accreted outside the effective radius 
and thus drive an episodical evolution of the half-mass radii, particularly below z = 1. In situ star 
formation and secular evolution processes contribute to 14^9^% and 16tn%, respectively. Overall, 
the simulated galaxies grow by a factor ^ 4 — 5 in mass and size since redshift z ^ 2. 

The short cooling time in the center of groups can foster a "hot accretion" mode. In one of the three 
simulated groups this leads to a dramatic rejuvenation of the central group galaxy at z < 1, affecting 
its morphology, kinematics and colors. This episode is eventually terminated by a group-group merger. 
Mergers also appear to be responsible for the suppression of cooling flows in the other two groups. 
Passive stellar evolution and minor galaxy mergers gradually restore the early-type character of the 
central galaxy in the cooling flow group on a timescale of ~ 1 — 2 Gyr. Although the average properties 
of central galaxies may be set by their halo masses, our simulations demonstrate that the interplay 
between halo mass assembly, galaxy merging and gas accretion has a substantial influence on the star 
formation histories and z = morphologies of central galaxies in galaxy groups. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: structure 



1. INTRODUCTION 

The study of galaxy groups is motivated by the fact 
that they host many galaxies in the local universe (Eke 
et al. 2004) and provide an environment in which galaxy 
interactions and merging are preferentially able to drive 
the morphological evolution of galaxies (e.g. Zabludoff 
& Mulchaey 1998). In addition, cluster galaxies likely 
experienced at some point in their history preprocessing 
in groups (e.g. Zabludoff 2002, Kautsch et al. 2008, Patel 
et al. 2009, Perez et al. 2009). Observations also indicate 
that bound, relaxed groups in the local Universe often 
host massive (> 10^^ M©) early- type galaxies at their 
centers (Zabludoff & Mulchaey 1998, Yang et al. 2009). 
The exploration of galaxy evolution in groups therefore 
provides the link between the study of isolated galaxies 
that preferentially populate lower mass halos (^ a few 
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10^^ ^0) and the analysis of member galaxies of the 
more massive galaxy clusters (> 10^^ ^o)- 

With the advent of large surveys it has now become fea- 
sible to identify and study the evolution in galaxy groups 
both at low (e.g. Weinmann et al. 2006a, Weinmann 
et al. 2006b, Yang et al. 2009, Guo et al. 2009, Carollo 
et al. 2009, Cibinel et al. 2009a, b,c) and high redshifts 
(e.g. Knobel et al. 2009) with an exquisite statistics. 
The study of early type galaxies in galactic halos has 
gained considerable attention (Meza et al. 2003; Naab 
et al. 2007, 2009), but the high computational demands 
necessary to resolve reliably the structural properties of 
individual galaxies pose a challenge for works that ad- 
dress the galaxy group scale (Kawata & Gibson 2003; 
Kobayashi 2005; Ohorbe et al. 2006, 2007; Kawata & 
Mulchaey 2008). Typically a spatial resolution of > 
1 kpc is reached (Kawata & Gibson 2003; Ohorbe et al. 
2006, 2007; Kawata & Mulchaey 2008) or the simulations 
are not fully cosmological (Kobayashi 2005). In this work 
we present a set of high-resolution simulations of galaxy 
groups that allow not only a study of the evolution of 
their central galaxies, but also enable us to investigate 
a rich satellite population within the virial radius of the 
group. Here, we put particular focus on the structural 
and kinematic properties of the central group galaxies 
and their most massive progenitors over time, relate their 
morphological properties to their assembly and gas accre- 
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tion histories, and analyze the evolution of masses, sizes 
and densities of the stellar component since z ^ 1.5. In 
particular, we address the following questions: 

1. Is the ACDM concordance model together with our 
current understanding of the baryonic physics able 
to reproduce massive galaxies that resemble those 
observed at the centers of z = galaxy groups? We 
perform our numerical experiments with the same 
code (Gasohne; Wadsley et al. 2004), with the same 
set of parameters, and at a resolution comparable 
to the one that has been used to study the forma- 
tion and evolution of disk and dwarf galaxies (Gov- 
ernato et al. 2007, Mayer et al. 2008, Governato 
et al. 2009a, Governato et al. 2009b). Therefore 
our simulations are not specifically tuned towards 
modeling massive galaxies in groups, but rather we 
apply a successful physical model to a yet unex- 
plored mass regime. 

2. What are the specific properties of the z ^ 2 pro- 
genitors of z = central group galaxies and how 
do they compare with the galaxies observed at those 
redshifts? One of the biggest challenges in obser- 
vational cosmology is to connect and relate galaxy 
populations at different times. Here the simple 
question is whether some of galaxies that we can 
observe at high redshifts, such as distant red galax- 
ies (Franx et al. 2003), restframe UV/optical se- 
lected galaxies (BM/BX Adelberger et al. 2004, 
Steidel et al. 2004), or BzK selected galaxies (Daddi 
et al. 2004), contribute to the precursor population 
of local central group galaxies, and, if so, what 
physical mechanisms play a role in the evolution 
of such systems from high redshift to the current 
epoch. 

3. Which physical processes contribute to the mass 
and size evolution of massive galaxies over the 
z ^ 2 to z = time span? Our simulations 
show that the progenitors of central group galaxies 
grow by a factor of 4-5 in mass and size between 
z ^ 2 and z = 0. Natural questions are there- 
fore (i) which physical mechanism, such as major 
and minor merging or star formation, dominates 
the mass and size evolution and (ii) how does the 
mass growth relate to the size evolution; in partic- 
ular, how steep is the size growth per unit mass? 
For example, minor merging has been suggested as 
the main driver in the evolution of passive, com- 
pact galaxies (e.g. Naab et al. 2007, 2009; Bezan- 
son et al. 2009) as it allows for a fast size growth 
per unit mass. Since we find that the z > 2 pro- 
genitors of the z = central group galaxies are star 
forming, their mass and size evolution may proceed 
differently. 

4. What is the impact of different assembly and gas 
accretion histories in groups at a fixed z = halo 
mass scale of ~ 10"^^ on the detailed physical prop- 
erties of their massive central galaxies? The as- 
sembly of the stellar component of the central 
galaxies is not simply a refiection of the assem- 
bly of the dark matter halo, but depends on the 
properties of the accreted baryons. As their halo 



TABLE 1 

COSMOLOGICAL PARAMETERS USED IN THIS WORK. 

Qrn ^A h ag n 

0.24 0.76 0.04185 0.73 0.77 0.96 

mass grows above 10^^ it is expected that these 
galaxies will shift from cold mode to hot mode ac- 
cretion (Keres et al. 2005; Dekel & Birnboim 2006; 
Brooks et al. 2009). However, the accretion mode 
may also be influenced by details of the environ- 
ment, such as by the properties of the surround- 
ing network of intergalactic gas fllaments (Keres 
et al. 2005). By z ~ all three groups have similar 
virial masses. They are therefore ideally suited for 
studying the impact of the rich variety of stochastic 
physical processes on the properties of their z = 
central galaxies. 

The paper is organized as follows. In sections 2 and 
3 we describe the simulation set-up and data analysis, 
respectively. In section 4 we present the properties of 
the simulated central galaxies over time and analyze the 
evolution of their masses, sizes and densities. We sum- 
marize our findings in section 5 and point out open ques- 
tions in section 6. In the appendix we describe how we 
deal with artificially enhanced central star formation (ap- 
pendix A), present our resolution test strategy (appendix 
C) and provide supplementary material (appendix D). 

2. SIMULATIONS 

The galaxy groups are selected from a DM-only simula- 
tion (Hahn et al. 2007) based on their Friends-of-Friends 
(FoF) masses (Geller & Huchra 1983; Davis et al. 1985) at 
redshift z = 1 using a standard linking length of 0.2. The 
employed FoF mass range is 0.8 - 1.2 x 10^^ M©. We note 
that while the original DM run uses WMAPl cosmolog- 
ical parameters (Spergel et al. 2003) our re-simulations 
are performed with a WMAP3 cosmology (see Table 1, 
Spergel et al. 2007) and therefore the mass selection cor- 
responds to a selection at redshift z ~ 0.7 in our cos- 
mology. Gl is an isolated group without any halo above 
1.5 X IO^'^MqH-^ within 5 Mpc/h at z = in WMAP3. 
The matter overdensity 6 = p/p— 1 within this radius 
is close to ((^ = 0.2). G3 on the other hand comprises 
a cluster of 1.6 x 10^^ MQh~^ and two massive groups of 
3.4 and 5.8 x IO^^MqH-^ within this radius {S = 15.0). 
Group G2 lies in-between the two extremes {6 = 1.4). 

The initial power spectrum has been generated with 
linger (Bertschinger 1995). Regions enclosing each 
galaxy group are refined with several layers of resolu- 
tion using graf ic-2 (Bertschinger 2001). The highest- 
resolution region is defined as the initial Lagrangian 
patch that contains the set of particles that enter at any 
point a sphere of radius R{z) around the particular group 
center. The group center is defined as the position of the 
particle with the highest density within the group halo 
(or its main progenitor). We use R{z) = 2 x Ryir{z = 0) 
fixed in comoving coordinates. The highest-resolution re- 
gion of each group is embedded into spherical Lagrangian 
patches of increasing radius and decreasing resolution. 
Gas particles are added only in the highest-resolution 
regions to reduce the computational costs. 

Our simulation strategy is as follows: We simulate 
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three galaxy groups (Gl, G2, G3) with SPH at "inter- 
mediate resolution" (see Table 2) down to z = 0. In 
addition, we simulate our fiducial group G2 at "high res- 
olution", i.e. 8 times higher mass and 2 times higher 
force resolution, down to z ~ 0.1 and denote this simula- 
tion as G2 — HR. We have resimulated this group also at 
a varying number of coarser mass and/or force resolution 
in order to study the impact on our results. We keep the 
softening length fixed in physical coordinates from red- 
shift to a high redshift, here 10, and fixed in comoving 
coordinates for higher redshifts, as found beneficial in 
simulations of disk galaxies (e.g. Governato et al. 2004, 
2007) or galaxy clusters (e.g. Borgani et al. 2006). We 
discuss the results of our resolution tests in appendix C. 

The simulations are performed with the TreeSPH code 
Gasoline (Wadsley et al. 2004) which is based on the 
parallel, multiple time stepping N-body code PKDGrav 
(Stadel 2001). The gravity tree opening angle is 0.525 
above redshift 2 and 0.7 below it. 

We use a standard radiative cooling scheme for primor- 
dial (metal-free) gas (Wadsley et al. 2004; Stinson et al. 
2006). We discuss changes expected due to metal cooling 
in appendix B. The simulations include a spatially uni- 
form UV background field (Haardt & Madau 1996). The 
UV background reduces the formation of galaxies with 
low masses (< lO^M©) and also implies a lower limit on 
the halo mass of the order of ~ 10^^ Mq above which 
accreted gas is shock heated (Brooks et al. 2009). 

We model the star formation and feedback following 
Stinson et al. (2006). Stars are formed in a probabilistic 
fashion if the gas density is larger than 0.1 cm~^, the 
gas temperature is lower than 15 '000 K, the gas is in an 
overdense {6 > 55) region and part of a convergent flow. 
The star formation efficiency is set to 0.05. The simula- 
tions incorporate supernova type la and type II feedback 
and mass loss of star particles due to stellar winds. Each 
star particle is treated as a single stellar population with 
Scalo IMF (Miller & Scalo 1979) and the parametrization 
of Raiteri et al. (1996) of the stellar tracks. The mod- 
elling of the type II supernovae makes use of the analytic 
blastwave scenario of McKee & Ostriker (1977). In par- 
ticular, a thermal feedback of 4 x 10^^ erg per supernova is 
injected into neighbouring gas particles which have their 
cooling shut-off for the time corresponding to the end 
of the snowplow phase of the blastwave. The cooling is 
not disabled for supernovae of type la. Over its lifetime 
subgrid stellar winds return a substantial fraction of the 
mass of a star particle (^ 40%) to its surrounding gas 
particles. 

3. METHODOLOGY 

We identify halos and galaxies and determine virial 
radii and masses with the help of the AMIGA halo finder 
(AHF, Gin et al. 2004; Knollmann & Knebe 2009). Halos 
extend out to the virial radius, which we define as the 
radius enclosing a mean matter density of 200 /Qrn times 
the mean density of the universe, or (for subhalos) out 
to the tidal radius, whichever is smaller. Galaxy masses 
are then estimated both in real space and in projection. 
We identify and extract the central galaxy of our group 
using the star particle with the lowest potential in the 
halo as center. We assign the central galaxy a stellar 
mass by measuring the mass of each stellar particle in a 
sphere of 20 physical kpc around the center. The choice 



of this radius is to some extent arbitrary. However, the so 
determined stellar masses change typically by less than 
10% when the radius is varied within a factor of 2. We 
show in Appendix C that these masses are well converged 
(up to a ~ 20%) below z = 2. We observe some residual, 
artifical, star formation even at z = and correct our 
masses for this effect (see appendix A). This lowers the 
stellar masses at z = by ~ 20% and by less at higher 
redshifts. 

Much of our analysis is done on mock images. The 
image projections are chosen to be along the axes of the 
moment-of-inertia tensor of the stellar mass within 10 
kpc. The (reduced) moment-of-inertia tensor T is de- 
fined as 

k 

where x\ is the i—ih. component of the vector pointing 
from the halo center to the /c-th particle position, is the 
distance from the halo center and rrik denotes the mass 
of the k-ih particle. The mock images have a pixelscale 
of 0.25 physical kpc per pixel, i.e. comparable to the 
gravitational softening of 0.3 kpc in our high resolution 
simulations G2 — HR. The mass (or light) of a given 
stellar particle is assigned to its nearest grid point (NGP) 
in the projected image. We use the NGP approach for all 
mock images except Fig.l, where we employ instead an 
adaptive cloud-in-cell method. In this case the mass (or 
light) of each particle is spread into a cube of a size that is 
proportional to the local interparticle distance. In either 
way, contours lines of constant surface mass density are 
created with SAOImage DS9. 

We define as effective radius the radius that includes 
half the stellar mass or flux within 20 kpc around the cen- 
ter of the respective galaxy. The estimate may include a 
correction according to the minimal star formation cor- 
rection scheme, see section A. For the z = central 
galaxies we additionally derive effective radii from Ser- 
sic or deVaucouleurs fits to the surface mass-density or 
surface brightness using the intensity values determined 
by running the IRAF task ellipse on mock images. The 
images are created by projecting along the minor axis of 
the moment-of-inertia tensor of the stellar component. 
In the employed weighted least square fit we only use in- 
tensity values for circularized radii between 3 x ebar and 
20 kpc and consider integer Sersic indices in the range 
1-8. We adopt the Sersic index with the lowest value 
as our best fit solution. We find similar results whether 
we use mass- or /-band flux-weighted effective radii and 
refer to the former unless mentioned otherwise. 

The star formation rates (SFR) that we quote refer 
to the stellar mass formed within 20 kpc around the 
main progenitor of the central group galaxy (excluding 
contribution from substructure) within a time-span of 
~ 100 Myr. Note that an "archaeological" measurement 
at z = would give much larger SFRs at earlier times 
because it measures the formation of all stars that are 
within the galaxy by z = 0, rather than the SFR of the 
most massive progenitor that existed at a given time. 

We calculate the luminosity of each stellar particle 
from its age and metallicity assuming a single stellar pop- 
ulation with Chabrier IMF (Bruzual & Chariot 2003). 
We measure the ^'^{g — r) and ^-^(r — i) colors and the 
^•^Mi absolute magnitude of the central galaxies at z = 
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TABLE 2 

Particle masses and spline softenings of the different resolutions. 









mstar 




^bar 


Resolution 


(Mo /i-i) 


(Mo h-^) 


(Mo h-^) 


(kpc /i-i) 


(kpc /i-i) 


intermediate 


3.7 X 10'^ 


7.9 X 10^ 


2.3 X 10^ 


0.73^ 


0.44^ 


high 


4.7 X 10^ 


9.9 X 10^ 


2.9 X 10^ 


0.36 


0.22 



Note. — The second column denotes the mass of the dark matter particles in the zoom-in region. The following 
columns indicate the initial mass of gas particles, the initial mass of spawned star particles, the gravitational softening 
of dark matter particles and the gravitational softening of baryonic (gas and star) particles, respectively. 



^For G2 a shghtly larger gravitational softening of 0.88 kpc 
has been used. 

^For G2 a slightly larger gravitational softening of 0.53 kpc 
has been used. 

in the z = 0.1 redshifted SDSS filter system of Blanton 
et al. (2003b) and the rest-frame SDSS r-band magnitude 
Mr. We follow Blanton et al. (2003b) in determining the 
colors and magnitudes from within a circular aperture 
with a diameter of 4 Petrosian radii. To avoid biases 
towards bluer colors from artificial star formation in the 
unresolved, central region of a 2: = central galaxy we 
fix the mass-to-light ratio within 2 x ebar to its value at 
2 X ebar- In addition, we estimate the amount of mass 
deposition due to artifical star formation in the central 
region and subtract a corresponding, excessive flux from 
our mock images within a circle of 2 ebar around the cen- 
ter. We note that the colors of the z = central galaxies 
do not include a correction for dust absorption. In a 
similar way we determine the colors and magnitudes of 
the main progenitors at higher redshifts (see Table 5 for 
the set of employed filter transmission curves). We clas- 
sify the most massive progenitors with the help of the 
fluxes measured within a projected radius of 8 physical 
kpc (~ f ), excluding the light from resolved satellites. 
We do not apply a mass-to-light ratio fix or a subtrac- 
tion of central flux, mainly because we expect central 
star formation to be a physical reality at high redshifts 
and the impact of additional artificial star formation is 
consequently small. Nonetheless, the quoted color and 
magnitude errors include variations arising from different 
viewing angles and from the difference between correct- 
ing or not-correcting for artificial star formation. Ab- 
sorption by dust may affect the colors and magnitudes 
of the gas-rich central galaxies at high redshifts and we 
therefore explore extinction corrections (Calzetti et al. 
2000) with Ay ranging from to > 1. At high redshifts 
{z > 2) we use as default value Ay = 0.8 which is the 
mean value found in the (almost) mass-limited sample 
of (Kriek et al. 2008a, b) and comparable to the median 
value in the SINS sample (Forster Schreiber et al. 2009, 
see also Genzel et al. 2006). Colors and magnitudes are 
given in the AB system, unless explicitly noted otherwise. 

In order to measure stellar line-of-sight rotation veloc- 
ity and velocity dispersion we first project the galaxies 
along the intermediate axis of their moment-of-inertia 
tensor. Then, we put a slit of extent and width of 20 kpc 
and 2 kpc, respectively, along the major axis of the pro- 
jected image and measure velocity and dispersion in 24 
bins along the slit. We proceed in an analogous way for 
gas velocities and velocity dispersions. The velocity dis- 
persion within the effective radius is then estimated from 
the mass- weighted average of all bins within 2 x ebar and 
3.5 kpc (~ i^eff) of the velocity dispersion profile. Sta- 



tistical errors are estimated from a bootstrapping anal- 
ysis. The statistical errors are typically much smaller 
than systematic effects that arise from using either ve- 
locity moments or Gauss-Hermite polynomials to fit ro- 
tation velocities and dispersion. Due to the latter we 
estimate that our velocities and dispersions are accurate 
to - 10%. 

Bulge-to-disk ratios of the central galaxies are mea- 
sured on mock images with GALFIT (Peng et al. 2002). 
We compute noise ( "sigma" ) maps based on the particle- 
per-pixel number in mock images assuming a Poisson 
statistics. We then run GALFIT on each central galaxy 
16 times covering the following cases: (i) edge-on vs. 
face-on projection, (ii) projected mass density vs. I-band 
flux density, (iii) a single component (Sersic profile) vs. 
a double component (Sersic profile + exponential disk 
profile) model, (iv) including or excluding the central 2 
softening lengths from the fit. While GALFIT very often 
produces good (x^ ~ 1) and sensible fits, in some cases it 
can be stuck in a non-optimal solution (either not a good 
fit (x^ > 2) or a clearly unphysical solution). In these 
cases we either changed our initial parameters, guided 
by visual inspection, or introduced some constraints that 
cut away the unphysical solution space. 

We compute the contribution to the growth of the stel- 
lar envelope by monitoring the stellar mass flux into the 
region between the effective radius and 20 kpc (our de- 
fined size of the central galaxy) among pairs of successive 
snapshots. The merger flux is the net flux of stars into 
this region that have formed at an earlier time outside the 
central galaxy. The stellar transfer flux is defined as the 
net flux of stars that formed in situ earlier on and are now 
crossing the mean effective radius between the two snap- 
shots. Stellar particles that form from on snapshot to the 
other in the shell enclosed by the effective radius and 20 
kpc define the in situ star formation flux. These fluxes 
represent smoothed quantities averaged over the time in- 
terval between successive snapshots. We checked that 
our results are not affected by our sampling frequency of 
the snapshots. An alternative methode to obtain fluxes is 
to apply the continuity equation using positions and ve- 
locities from a single snapshot. However, unless the time 
resolution is exquisitely high this approach is not well 
suited to measure the average flux from non-continuous, 
clumpy accretion events and galaxy mergers for a given, 
fixed radius. 

In this paper we consider a merger between two galax- 
ies as major (minor) if the stellar mass ratio IZ is larger 
(smaller) than 1:3.5. Sometimes stars form in places that 
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are not identified as bound halos by our halo finder. In 
this case we say that these star particles originate from 
unresolved objects. In order to determine whether a par- 
ticular star particle is added to the central galaxy by 
a major merger, minor merger or from an unresolved 
source we proceed as follows. First we identify the last 
satellite halo to which this star particle belonged before it 
finally merged with the central galaxy. We keep track of 
whether the star particle is stripped from this last satel- 
lite halo before it finally enters a 20 kpc boundary around 
the central galaxy. If the star particle first appears in a 
resolved satellite located within a 20 kpc radius around 
the central galaxy, we say that the star particle forms by 
a merger-induced star formation event in the satellite. 
In case no such satellite halo can be identified we dis- 
tinguish between star particles that form within 20 kpc 
around the central galaxy (in situ star formation) and 
star particles that form outside this radius (unresolved 
origin). 

4. RESULTS 

In this section we address the four major questions 
that we raise in the introduction, namely: (i) How do 
predictions from ab-initio, cosmological, hydrodynamical 
simulations compare with observations of nearby central 
group galaxies? (ii) How do early-type central group 
galaxies at z = relate to galaxies observed at z ~ 2? 
(iii) How do the main progenitors of central group galax- 
ies evolve in their most fundamental parameters, namely 
mass and size? and (iv) How does the assembly history 
affect the evolution of central group galaxies? 

4.1. Overview of the groups and their central galaxies 

In Fig. 1 we show how the three groups Gl, G2 and 
G3 and their progenitors evolve from 2;~2.5toz = 0. 
At z > 2 the group progenitors are surrounded by fil- 
aments of cold gas in agreement with a cold accretion 
picture, e.g. Keres et al. (2005). These filaments disap- 
pear at z < 1.5. The (shock-heated) hot gas is confined 
to halos and filaments, i.e. regions of high dark matter 
density. The contours of constant surface mass density of 
the hot gas can have very irregular shapes at z > 0.6, but 
at lower redshift they become more regular. Stars and 
galaxies form in the places of high density in the cosmic 
filaments. A "central galaxy" forms at the bottom of the 
potential well of each group progenitor. Other galaxies 
approach and fall into the forming galaxy groups. Some 
of these galaxies merge quickly with the central galaxy, 
while others continue to orbit for a long time. At z ^ 
the virial masses of all three groups are in the range 
1 — 2 X IO^^Mq, although their amount of substructure 
varies substantially. The number of satellites with a stel- 
lar mass in excess of 10^° M© is 2 (Gl), 8 (G2), 8 (G3, 
ah at z = 0) and 13 (G2 - HR, at z = 0.13). 

Zooming in to galactic scale we show in Fig. 2 and 
Fig. 3 Bessel B, R and I band composite images of the 
central galaxies. Clearly, the gross evolution proceeds in 
all three cases in a similar fashion including a substan- 
tial change in color (from blue to red), morphology, size 
(from < 1 kpc to several kpc) and cold gas content (from 
gas-rich at z = 2.5 to gas-poor at z = 0). The most 
notable differences between the three groups are: (i) the 
galaxies show a varying degree of disky-ness at z = 0, 
(ii) at z ~ 0.6 group G2 harbors a massive spiral galaxy 



and a substantial reservoir of cold gas, (iii) group G2 
and its central galaxy appears more massive and evolved 
especially at z > 1.5. We now proceed with a more quan- 
titative analysis of the properties of the central galaxies. 
In the following we will, unless noted otherwise, include a 
minimal correction scheme for artificial, central star for- 
mation as described in appendix A. There we also discuss 
implications of more aggressive correction schemes. 

4.2. The properties of the central galaxies at z=0 

In this section we show that the central objects resem- 
ble early-type galaxies based on various criteria. We also 
compare the properties of the simulated central group 
galaxies with observations of central group galaxies or of 
massive early-type galaxies in general, and with results 
from previous simulations. 

Galaxies can be split into early- and late-types accord- 
ing to different criteria, e.g. based on colors, Sersic index. 
Bulge-to-disk ratios, rotational support, or morphology. 
In addition, indirect indicators can be used - such as 
masses (the most massive galaxies are early-type galax- 
ies) or gas fractions (early-type galaxies are often gas- 
poor). Using various criteria we now show the early- type 
nature of the simulated central group galaxies. 

The red ^'^{g — r) and ^-^(r — i) colors of the central 
group galaxies, see Table 4, are indicative of an aging 
stellar population. Figure 4 shows that surface mass den- 
sity profiles of the central galaxies are close to a deVau- 
couleurs profile. When fit with a Sersic law, the indices 
are 4, 5 and > 8, respectively. Analyzing the mass den- 
sity or I-band images with GALFIT (Peng et al. 2002) 
we find that Gl and G3 are very well fitted (x^ ~ 1 — 1.3) 
by a single n > 3 Sersic profile (independent of inclina- 
tion and of whether or not we mask the central region). 
G2 can also be fitted with a single high-n Sersic profile 
as long as the central two softening lengths are excluded. 
Alternatively, it can be well fitted with a two-component 
model in which either (a) a small disk (scale length ~ 0.7 
kpc) is embedded in a more extended spheroid (effective 
radius of 3.4 - 5.5 kpc depending on whether or not the 
central region is masked), or (b) a central spheroid (1.6- 
1.8 kpc) is embedded in an extended, but faint disk (scale 
length ~ 6 — 7 kpc). The first model is favored in edge-on 
projections, the second in face-on projection. The bulge- 
to-disk ratios of Gl and G3 are naturally very high and 
are typically > 3. For G2 we obtain a bulge-to-disk ratio 
of 3 (edge-on) and 1.6 (face-on). Also a visual inspec- 
tion leads to the conclusion that the central galaxies are 
of early-type. The central galaxies in the groups Gl and 
G3 have an elliptical morphology while the central galaxy 
of G2 clearly shows some disk component. Spiral arms 
in the disk are revealed in its high resolution equivalent 
G2-HR at 2: ~ 0.6, but disappear later on. At z ~ 
the central galaxy of group G2 resembles a SO or E/SO 
galaxy, see Fig. 2 and Fig. 3. 

In Figure 5 we plot the stellar rotation velocity and 
velocity dispersion along the intermediate axis of the 
moment-of-inertia tensor of three groups Gl — G3. Pro- 
jecting along this axis results in the largest line-of-sight 
rotation velocity and we confirm that the stellar angu- 
lar momentum is well aligned with the minor axis of the 
stellar component of the galaxy. The figure demonstrates 
that the central galaxies are supported by velocity dis- 
persion. We note that the central galaxies in the groups 
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Fig. 1. — The appearance of the simulated galaxy groups as function of redshift. The size of each image is 1.2 Mpc x 1.2 Mpc. The four 
columns correspond to z = 2.5, z = 1.5, z = 0.6 and z = (z = 0.1 for G2 — HR), while the rows correspond to the different groups CI, 
G2, G3, and G2 — HR (from top to bottom). Color coded are dark matter in blue (from 3.6 to 1460 Mq pc~^), cold gas (here the gas 
with T < 2.5 X 10^ K) in green, stellar matter in yellow (both from 0.9 to 365 Mq pc"^), and hot gas (here the gas with T > 2.5 x 10^ 
K) as red surface mass isocontours (3 contour lines at 1, 4.5 and 20 Mq pc~^). The white circle shows the virial radius of each group at 
the indicated redshift. 



G2 and G3 {v/a ^ 0.7 — 0.9) show a significant rotation 
component, while the central galaxy in Gl is practically 
non-rotating. Interestingly, the last major mergers that 
the central galaxies in the groups Gl and G3 undergo 
are between gas-poor, velocity dispersion supported stel- 
lar systems. We explain the relatively high rotational 
support of the remnant in the group G3 by the special 
orbital properties of merging galaxies. Indeed, a kine- 
matic re-analysis of the binary merger simulations that 
have been presented in Feldmann et al. (2008) shows that 
a gas-poor major merger on an eccentric orbit, e.g. with 
an apo-to-pericenter ratio of 6:1, typical for orbits of DM- 
halos in bound environment such as clusters and massive 
groups (Ghigna et al. 1998), can also lead to v/a of or- 
der unity. Note that typical binary merger experiments 
assume parabolic orbits which are expected for merg- 



ers between field galaxies in isolated halos (Khochfar & 
Burkert 2006), and which convert a lower amount of an- 
gular momentum into the spin of the merger remnant 
(Cox et al. 2006). 

The high masses of the simulated galaxies (~ 4 x 
10^ within a sphere of 20 kpc radius) and the 
scarcity of cold gas (< 10^ Mq of gas colder than 3.2 x 10^ 
K within a sphere of 20 kpc radius) are also consistent 
with the expectation of being early-type galaxies. We 
summarize the properties of the central galaxies at z = 
in tables 3 and 4. We conclude that based on visual ap- 
pearance, lack of cold gas and rotational support, color, 
surface-mass profile and bulge-to-disk ratio, the central 
group galaxies resemble early-type galaxies. 

Our simulated central galaxies have larger masses and 
velocity dispersions and/or smaller sizes (Meza et al. 
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Fig. 2. — Mock images of the central group galaxies (and their main progenitors) as function of redshift. Each image is 20 physical kpc 
across and shows a face-on projection corresponding to the minor axis of the moment-of-inertia tensor of the stellar mass within a sphere of 
10 kpc radius. The four columns correspond to 2; = 2.5, z — 1.5, z = 0.6 and z = Q {z = 0.1 for G2 — HR). The first three rows correspond 
to Gl, G2, and G3 (from top to bottom) while the last row shows the appearance of the central galaxy in the high resolution re-simulation 
G2 — HR. When assessing the images note that (i) the projections are along the line-of-sight of the virial radius of the galaxy group, (ii) 
the image orientations are arbitrary, hence the displayed disk orientations are not directly comparable, (iii) the addition of high frequency 
modes in the high resolution re-simulation and low-level numerical noise can lead to subtle timing and positioning differences at small scales 
and (iv) the central galaxy in group G2 (G2 — HR) undergoes a substantial morphological evolution between z = 0.2 and z = 0. The RGB 
color channels of the images correspond to the surface brightness in the restframe Bessel B, R and I filterbands, respectively, and range 
from 13.5 mag arcsec"^ to 22 mag arcsec"^. Green contours indicate column densities of cold gas (T < 3.2 x lO^K) corresponding to 1, 
10, 100, and 1000 Mq pc~^. The green circle has a radius of 2 gravitational softenings Cbar? indicating the resolution limit. The white 
circle encloses the central 2 kpc. 



2003), when compared with galaxies observed in the 
nearby Universe (Shen et al. 2003, van der Wei et al. 
2008, Fig 3). The compactness of our objects is less ex- 
treme, however, than e.g. in Meza et al. (2003) (the effec- 
tive densities being an order of magnitude lower), which 
we ascribe to our efficient supernova blastwave model 
and, partially, to the increased resolution. In order to 
follow the average mass-size relation of local early-type 
galaxies (Shen et al. 2003) or local central group galaxies 
(Guo et al. 2009) the simulated central galaxies would 
need to reduce their mass by a factor of 2-3. This is also 



the amount necessary to put the galaxies right onto the 
red-sequence in the ^'^{g — r) and ^-^(r — i) vs. ^-^Mi 
color-magnitude diagrams of nearby galaxies (Blanton 
et al. 2003a). Similar to what is obtained from simu- 
lations of central galaxies in galactic halos (Naab et al. 
2009) our objects appear to be too massive, by about a 
factor of 2-3, for their halo mass when compared with 
weak lensing observations (Mandelbaum et al. 2006), or 
compared to SDSS groups (Yang et al. 2007, 2008; Guo 
et al. 2009), although this comparison is hampered by 
the fact that masses are determined in observations and 
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Fig. 3. — As Fig. 2 but showing instead an edge-on projection corresponding to the intermediate axis of the moment-of-inertia tensor of 
the stellar mass. 



TABLE 3 

Structural and kinematic properties of the central group galaxies at z = 0. 

Group (1012 Mo) (kpc) (IQIOMo) (kpc) (kpc) (kpc) n^^rsic (km/s) Vrot/(Tefi 



Gl 11.0 447 39.6 ±3.6 4.0 ± 0.5 2.75±0.1 2.50±0.1 5 299±20 0.20±0.05 

G2 12.0 459 44.5 ±3.6 3.7 ±0.5 3.81±0.2 3.81±0.2 4 326±20 0.79±0.05 

G3 15.6 502 42.9 ±3.6 3.2 ± 0.4 2.57±0.1 1.57±0.1 > 8 323±20 0.86±0.06 

Note. — The first column lists the group name. Further columns denote: the virial mass and virial radius of the group, the stellar 
mass M* within 20 physical kpc, the three-dimensional radius that contains half of M*, the half-mass radius from a deVaucouleurs-fit, the 
half-mass radius from a Sersic fit, the Sersic index, the stellar velocity dispersion within the effective radius, the ratio of stellar rotation 
velocity and velocity dispersion. Errors in M* and i^eff ^ire derived from varying the applied star formation correction in the range 6-16 
Mq/jt. Errors in R^^^^^'^ and R'^^^^ are formal fit errors. Errors in the kinematic properties are derived from a bootstrapping error 
analysis of the line-of-sight velocity data. Masses and sizes derived without the minimal star formation correction scheme can be found in 
Table Al. In Table CI we compare the results of the simulations G2 and G2 — HR at z = 0.1. 
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TABLE 4 

Photometric properties of the central group galaxies at z = 0. 



etrosian DPetrosian 
^eff 

Group (kpc) (kpc) ^-^(g-r) ^-^(r-i) ^-^Mi Mr 



Gl 6.9 ±1.3 3.0 ±0.5 0.87 ± 0.02 0.36 ± 0.02 -23.51 ±0.11 -23.34 ±0.11 
G2 6.7 ±1.8 2.6 ±0.6 0.81 ± 0.02 0.34 ± 0.02 -23.87 ±0.10 -23.72 ±0.10 
G3 4.8 ±0.5 2.1 ±0.4 0.86 ± 0.02 0.35 ± 0.02 -23.54 ±0.11 -23.38 ±0.11 

Note. — The first column lists the group name. Further columns denote: the Petrosian radius in the SDSS r-band, the radius containing 
half of the light within an annulus of 2 Petrosian radii, ^'^(g — r) color, ^ ^(r — i) color, and ^-^i band magnitude in the redshifted SDSS 
u,g,r,i,z system (Blanton et al. 2003a) and the restframe SDSS r-band magnitude. Errors include variations due to different viewing angles 
and due to a star formation correction in the range 6-16 M0/yr. 




Fig. 4. — (From left to right) The surface mass density at z = of the central objects in the group Gl, G2 and G3. The best-fit Sersic 
profile (blue line) with n G {1, ... ,8}, the best-fit deVaucouleurs profile (magenta line) and the derived effective radii are indicated. Fit 
solutions inside (outside) the fitting range are shown as solid (dashed) lines. The shaded region, which corresponds to the (projected) 
central region within 3 gravitational softening lengths Cbar? has been excluded from the fit. 
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Fig. 5. — The stellar rotation velocity (red) and velocity dispersion (black) along the axis which maximizes rotation velocity of the central 
objects in the three groups Gl, G2 and G3 at 2; = (from left to right). The dashed curves show the velocities and velocities dispersions 
using moments of velocity, while the solid curves are calculated from fitting a 5th order Gauss-Hermite polynomial to the line-of-sight 
velocity distribution. The dark shaded area is affected by numerical resolution and excluded in the analysis. We quote in table 3 the 
velocity dispersion measured within the light shaded area. 
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simulations in a very different manner and a systematic 
(but not statistical) uncertainty of about a factor of 2 is 
not unlikely. Finally, in Fig. Dl we plot the integrated 
mass function of Baldry et al. (2004). If we assume a 
1:1 relationship between central galaxy and their host- 
ing dark matter halo we can associate the number densi- 
ties of galaxies with the number densities of their parent 
dark matter halos. The central galaxy in a ~ 10^^ Mq 
dark matter halo Mo & White (2002) is then expected 
to have a stellar mass of ~ 2 x 10^^ M©, i.e. a factor of 
2 smaller than the masses of the central group galaxies 
in the simulations. Taken together these results indicate 
that the masses and sizes of the central group galax- 
ies in high-resolution cosmological simulations are still 
somewhat biased towards larger masses and/or smaller 
sizes, reminiscent of the size problem in simulations of 
disk galaxies. This is not entirely unexpected. The size 
problem is usually associated with the loss of angular 
momentum due to inefficient resolution at high redshift 
when the low-mass progenitor galaxies are resolved by 
only a small number of particles and with the particular 
feedback implementation (see Mayer et al. 2008 for a re- 
cent review). Hence, it is likely that our simulations are 
also affected to some extent by this problem. Nonethe- 
less, it should be stressed that reducing the stellar mass 
by a factor 2-3 would bring our galaxies in good agree- 
ment with observations and we note that this factor is 
likely not much larger than the systematical uncertainties 
in the masses determined from observations. To summa- 
rize: Our simulations produce massive early-type galax- 
ies at the group centers that compare reasonable well 
with observations although they are somewhat biased 
w.r.t. masses, sizes, colors and/or magnitudes. Possibly 
a combination of higher resolution and an increased star 
formation threshold compatible with the cold molecular 
phase and/or additional feedback mechanisms such as 
AGN feedback could resolve these remaining discrepan- 
cies (Kravtsov 2003, Robertson & Kravtsov 2008, Saitoh 
et al. 2008, Governato et al. 2009b). 

4.3. The z ^ 2 progenitors of z = central group 
galaxies 

The main progenitors at 2; > 2 of our studied central 
group galaxies have very compact (< 1 kpc) stellar com- 
ponents and given our spatial resolution we cannot re- 
liably estimate their structural or kinematic properties. 
However, global properties such as total mass or lumi- 
nosity can be robustly measured. We therefore classify 
the main progenitors at high redshift s based on colors, 
SFR and cold gas content. In addition, we compare their 
halo masses with observational estimates and clustering 
measurements. 

Frequently used classification and selection schemes for 
high redshift galaxies are the BzK color-color classifica- 
tion (Daddi et al. 2004), the BM/BX selection (Adel- 
berger et al. 2004; Steidel et al. 2004) and the selection as 
distant red galaxies (DRG) (Franx et al. 2003). In Fig. 6 
we show the evolution of the colors of the main progen- 
itors in the BzK plane between z = 3 and z = 0.5. At 
z > 2 the galaxies fall into the star forming BzK (sBzK) 
regime consistent with their high star formation rates of 
20 — 6OM0 yr~^ and specific star formation rates of > 0.2 
Gyr~^. The colors atz~2 (5 — z^l,z — Kg ~ 2) are 
typical of star forming galaxies observed at those red- 



shifts, e.g. Forster Schreiber et al. (2009). At ^ - 1.3 
the galaxies cross the line in the BzK plane that is typi- 
cally used to select either z < 1.4 galaxies or z > 1.4 star 
forming galaxies. None of our galaxies classifies at any 
time as a proper passively evolving BzK (pBzK) galaxy, 
although the central galaxies in the groups Gl and G2 
come temporarily close to the pBzK area at z ~ 1.5. 
Kriek et al. (2008a) have shown that massive i^-bright 
galaxies observed at z = 2 — 3 divide into (i) red, post- 
burst galaxies that form a red sequence in a mass vs. 
{U - 5)rest-frame,Vega plauc and (ii) bluc star forming 
galaxies with small star formation timescales. The blue 

{U - 5)rest-frame,Vega ColorS of OUr objcctS at Z > 2 (- 

0.2 to 0) would put them into the second category. At 
z ^ 1.5 — 2^ however, the central galaxies of the groups 
Gl and G2 are redder (> 0.1) and could be classified 
as post-burst galaxies. In Fig. D2 and Fig. D3 we show 
how the progenitor galaxies would classify according to 
the DRG and BM/BX scheme, respectively. These se- 
lection schemes are somewhat more susceptible to dust 
extinction and we find that our classification outcome 
depends on the adopted extinction value. For Ay < 1 
(and in particular for the assumed value Ay = 0.8, see 
section 3) our simulated progenitors classify as BX galax- 
ies in the redshift range z ^ 2.7 and z ^ 2. The less 
massive progenitor galaxies (the progenitors of the cen- 
tral galaxies in the groups Gl and G3) qualify also as 
BM galaxies below z = 2 down to z ^ 1.7 (Gl) and 
z ^ 1.5 (G3). Only if the progenitors were strongly dust 
obscured {Ay > 1.3) they would qualify at z > 2.3 as 
distant red galaxies (Franx et al. 2003). We summarize 
the photometric properties in Table 5. 

Masses and star formation rates of the central galaxies 
at z ~ 2 are summarized in Table 6. Consistent with 
the blue colors are the significant star formation rates, 
see Fig. 12, which are rapidly declining with time. At 
z ^ 2 they still amount to 20-60 Mq yr~^. The progen- 
itor galaxies also host a significant reservoir of cold gas 
(5x 10^ Mq at z ~ 2 — 2.5 within a 10 kpc radius) which is 
typically arranged in either a gas disk of ~ 3 kpc radius 
(in G2 and G3) or it has a more irregular morphology 
(Gl). When we look at the kinematic properties of this 
cold gas disk in our best resolved object (G2 — HR) we 
measure a line-of-sight velocity dispersion of 200 km/s 
and a two times higher rotation velocity in the gas disk. 
We note that compared to e.g. the star forming galaxy 
reported in Genzel et al. (2006) the cold gas masses in our 
objects are significantly lower (by about an order of mag- 
nitude) and our disk rotation speed is larger (factor of 2) 
and reaches its maximum value at smaller radii. The lat- 
ter is likely a consequence of the fact that the simulated 
galaxies are rather compact with stellar half-light radii 
below 1 kpc. We note that also about 1/3 of massive, star 
forming galaxies seen at 1.5 — 2.5 are compact and 
have high velocity dispersions (Forster Schreiber et al. 
2009). On the other hand, the size difference at z ^ 1.5 
between the simulations G2 and G2 — HR (see Fig. 8) 
indicates that our galaxies might suffer to some degree 
from artificially enhanced angular momentum loss, often 
seen in simulations of disk galaxies, despite the fact that 
our model invokes relatively energetic supernova feed- 
back. Clearly, this issue needs further studies at higher 
resolution. Overall, we find that the colors and star for- 
mation rates of the simulated galaxies match the proper- 
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(B-z)ab 

Fig. 6. — The BzK colors of the main progenitors of the central 
group galaxies within a projected radius of 8 kpc between z = 3 
and z = 0.5. The different lines corresponds to Gl (green), G2 
(red), G3 (blue) and G2 — HR (magenta). The gray symbols indi- 
cate the change of the colors due to different extinction corrections 
(Calzetti et al. 2000) with Ay ranging from to 1.3. The color 
evolution is shown for the default value Ay =0.8. Filled (empty) 
symbols denote z > 1.4 (z < 1.4). Squares, triangles and cir- 
cles indicate that the specific star formation rate within 20 kpc 
is > 0.5 Gyr-i, between 0.2 and 0.5 Gyr^^ or below 0.2 Gyr^^, 
respectively. The black star indicates the BzK colors of a typical 
z 2 star forming galaxy (Genzel et al. 2006). The error bar at 
the top shows the maximum changes that result from changing the 
projection direction and when the minimal star formation correc- 
tion scheme is applied. Satellite galaxies that happen to lie along 
the line of sight are excluded because their presence can affect the 
overall colors. 

ties of optically/UV (BM/BX) selected z ^ 2 star form- 
ing galaxies (Genzel et al. 2006, 2008) or that of sBzK 
galaxies (Daddi et al. 2004), but not that of high redshift 
"red-and-dead" galaxies. 

Measurements of the angular correlation functions al- 
low to determine the typical halo masses in which sBzK 
galaxies of a given magnitude reside (Kong et al. 2006; 
Hayashi et al. 2007; Blanc et al. 2008; Hartley et al. 
2008). The parent halo mass increases rapidely with K- 
magnitude Hayashi et al. (2007): by more then 2 orders 
of magnitude per 2.5 mag in K brightness. The brightest 
sBzK galaxies {Ks < 21) thus populate halos compa- 
rable to that of pBzK selected galaxies. Typical halo 
masses of sBzK galaxies reported in the literature are: 
2.8 X lO^^M© for a i^^ < 23.2 sample (Hayashi et al. 

2007) , 6 X lO^^Mo for a i^^ < 23 sample (Hartley et al. 

2008) and - lO^^M© for a Ks < 22 sample (Kong et al. 
2006; Hayashi et al. 2007; Blanc et al. 2008). The halo 
masses (0.5 — 1.8 x lO^^M©) and i^^-band magnitudes 
(21.5 — 22.4) of the simulated central galaxies are consis- 
tent with their identification as sBzK galaxies of inter- 
mediate brightness and star formation intensity z ^ 2. 

4:A. Evolution of masses, sizes and densities 

4.4.1. Total masses and effective radii 

The left panel of Fig. 8 presents the growth of stellar 
masses and effective radii of the central group galaxies as 



function of redshift. Mass and size evolution are closely 
linked and periods of slow/fast mass build-up correspond 
to periods of slow/fast size growth. The central galaxies 
in the groups Gl and G3 undergo at least two major 
galaxy mergers (z ~ 1, z ~ 0.4 in case of Gl and z ~ 0.8, 
z ^ OA in the case of G3). Two of the 4 major mergers 
occur between rotation-dominated, disky, non gas-poor 
(/gas/stellar of fcw percent) galaxies, while the other 
two mergers take place later between an already gas- 
poor central galaxy and another gas-poor, non-rotation 
supported companion. The central galaxy in G2 {G2 — 
HR) does not experience any major mergers below z = 4. 

The overall behaviour of the mass-size evolution is 
shown in the right panel of Fig. 8. It roughly follows 
a relation R oc M". The large jumps with a < 1 are 
major mergers, consistent with estimates from binary 
merger simulations (Boylan-Kolchin et al. 2006). Peri- 
ods of minor merging and star formations show a con- 
tinuos growth in mass and effective radius with a > 1. 
We perform a robust linear regression (with a bisquare 
weighting scheme) in the log (M) — log (i^eff ) plane in or- 
der to determine the average value of a over the redshift 
range z ~ — 1, see Table 7. We now summarize the 
result of this analysis. 

The two major merger below 2; < 0.5 have a = 0.25 and 
Q/ = 0.90. In phases without major merging activities 
we identify three mechanisms that drive significant size 
growth at small or only moderate mass growth. These 
processes are (i) minor merging, (ii) non-central star for- 
mation and (iii) a redistribution of either pre-accreted 
or pre-formed stellar material. The latter process may 
originate in a physical mechanism, such as tidal heating 
due to orbiting satellites, and/or is caused by spurious 
numerical effects. To assess whether the latter is the case 
we compare the intermediate resolution G2 and its high 
resolution analogue G2 — HR and find a resolution de- 
pendence of a. More precisely, at higher resolution the 
central galaxy seems to grow slower in size for a given 
mass {a ^ 1) compared to the corresponding simula- 
tion at intermediate resolution (a ~ 2). This behavior 
is partially explained by the fact that compared with G2 
the central galaxy in G2 — HR is larger at high redshift 
{z > 1) but of similar size at z ~ 0. It is clear that 
further work is necessary in order to tie down all resolu- 
tion dependent effects that potentially contribute to this 
difference. 

We conclude that major mergers alone typically result 
in a slow size growth (a < 1). Group G2 is subject to 
both substantial minor merging and star formation, and 
experiences a faster size growth (a ~ 1 in the highest 
resolution simulation). Finally, groups Gl and G3 show 
an even faster size evolution at nearly constant mass in- 
between major mergers. At least parts of this evolution 
could be of spurious numerical origin. However, we do 
not exclude that, e.g., heating by tidal shocks from or- 
biting satellites could contribute to this redistribution of 
stellar mass and thus drive a very fast size evolution, see 
section 4.4.3. 

In Fig. 7 we decompose the stellar mass found in cen- 
tral galaxies at given redshift into mass accreted by ma- 
jor/minor merging, mass produced by in situ star forma- 
tion in the central galaxy, mass produced in the satellite 
while it merges with the central galaxy, mass accreted 
after being stripped from a satellite and mass accreted 
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TABLE 5 

Photometric properties of the central galaxy progenitors at 2; ~ 2. 



Group 


z 




vjr — /V 


7^ 


1 K 

J s ~ J^s 


ID — Z 


Z — i\s 




Gl 


2.4 


0.57 


0.25 


24.3 


-1.7 


1.20 


0.58 


1.91 


22 




2.0 


0.40 


0.16 


23.1 


+1.1 
-1.6 


0.62 


0.37 


1.45 


21 ^+0-^ 
^-•-•^-0.6 




1.5 


0.86 


0.98 


24.5 


+0.9 
-1.4 


0.96 


2.52 


2.50 


20.8lg-J 


G2-HR 


2.4 


0.57 


0.22 


23.2 


-1.6 


1.02 


0.49 


1.55 


21.5ir. 




2.0 


0.86 


0.56 


24.9 


+1.0 
-1.7 


1.11 


1.50 


2.71 






1.5 


0.52 


0.87 


24.6 


+0.9 
-1.5 


1.10 


2.43 


2.69 


20.7ig-J 




2.4 


0.61 


0.26 


24.2 


-1.7 


0.99 


0.59 


1.6 


22 4+^-^ 




2.0 


0.51 


0.25 


23.6 


+1.1 
-1.6 


0.62 


0.53 


1.53 


21 9+°-^ 




1.5 


0.17 


0.22 


23.0 


+0.9 
-1.5 


0.77 


0.76 


1.81 


on 7+0.3 

' -0.5 



Note. — The first columns list the group name and the redshift of the simulation. The next columns contain the colors and magnitudes 
(all in the AB system) of the main progenitors of the z = central group galaxies in various filter bands used for selecting high-redshift 
galaxies. Un, G and 1Z filter bands are extracted from Fig. 1 of Steidel et al. (2003) and compare well with the transmission curves of 
(Adelberger et al. 2004). Jg is a slightly modified J filter (Labbe et al. 2003). The Jg and Kg, 2 transmission curves are downloaded from the 
ISAAC (VLT) website. We further use Bessel B, z (Subaru telescope), Kg^i (Kitt Peak 4-m telescope) filter bands. We use the Ks,i filter 
in order to compute z — K colors and the Ks,2 filter for calculating Js — Kg colors. Colors and Kg magnitudes are affected by much less 
than 0.1 mag when switching between Kg^i and Kg, 2 filter bands. Colors and magnitudes are measured in an observed frame and within a 
projected radius of 8 kpc (~ 1^^) roughly mimicking the observations of Franx et al. (2003); Labbe et al. (2003); Steidel et al. (2003, 2004); 
Hayashi et al. (2007). An extinction law (Calzetti et al. 2000) with Ay = 0.8 is assumed. The quoted error in the magnitude results from 
variations of Ay in the range — 1.3 which dominate the uncertainties in the magnitudes and colors. Since colors are strongly correlated 
with each other we do not quote corresponding uncertainties for them. Instead we refer to Fig. 6, Fig. D2 and Fig. D3 that demonstrate 
the impact of varying Ay on apparent colors and magnitudes. At fixed Ay variations in color and magnitude due to the chosen projection 
directions and the amount of SF correction in the minimal correction scheme (see Appendix C) typically amount to ^0.1 mag. 



TABLE 6 

Structural properties of the central galaxy progenitors at 2; ~ 2. 







Mvir 


-Rvir 


Mtot 


Mstar 


Mgas 


-^cgas 


SFR 


Group 


z 


(lOi^Mo) 


(kpc) 


(IOIOM0) 


(IOIOM0) 


(1O1OM0) 


(IOIOM0) 


(M0 yr-l) 


Gl 


2.4 


0.69 


53 


16.7 


4.9 


1.2 


0.39 


38 




2.0 


1.04 


68 


17.8 


5.7 


1.0 


0.35 


26 




1.5 


1.38 


89 


23.2 


9.9 


0.6 


0.19 


16 


G2- HR 


2.4 


1.52 


69 


27.1 


8.6 


1.4 


0.41 


48 




2.0 


1.82 


82 


28.2 


9.5 


0.8 


0.18 


22 




1.5 


2.55 


110 


31.9 


13.0 


0.6 


0.07 


14 


G3 


2.4 


0.49 


47 


11.8 


2.1 


1.3 


0.62 


14 




2.0 


0.61 


57 


16.9 


6.2 


1.6 


0.71 


60 




1.5 


1.28 


87 


21.0 


8.5 


1.1 


0.46 


26 



Note. — The first columns list the group name, redshift, virial mass and physical virial radius of the parent halo of the main progenitors 
of the z = central group galaxies. The next columns show total mass, stellar mass, gas and cold gas (T < 3.4 x 10^ K) mass within a 
sphere of 10 kpc radius (~ 1^^) around the object's center. The last column contains the star formation rate within a spherical radius of 20 
kpc around the progenitor galaxy, but excluding contributions from substructure. Cited stellar masses and SFR are corrected according to 
the minimal correction scheme which changes the stellar mass by < 10%. The mass subtracted from the stellar component is not added to 
the gas masses. The effective radius of the stellar component is not resolved, but smaller than ~ 1 kpc. 



that has been formed in an unresolved substructure out- 
side the central object. As the left panel shows only a 
small amount of stellar mass is accreted smoothly, e.g. 
either after being stripped from a satellite or after be- 
ing formed in an unresolved substructure. Most of the 
stellar mass is accreted by major or minor mergers or 
is produced by in situ star formation within the central 
galaxy. Group G2 which does not experience a single ma- 
jor galaxy merger below z = 4 forms 70—90% of its stellar 
mass in situ, while both Gl and G3 accrete about 70% of 
their final mass by merging. Minor merging contributes 
at the 5 — 25% level to the mass build-up of the central 
group galaxies. The right panel shows that the contribu- 
tion of in situ star formation is overall reduced compared 
to the effect of merging for 2; < 1, indicating that the star 
formation rate drops faster than the merger rate, as pre- 
dicted by semi-analytic modeling (Guo & White 2008). 
It is clear that a larger sample of galaxies with represen- 



tative merger and gas accretion histories will be required 
in order to pin down in a statistical manner the contribu- 
tions from the various processes (major merging, minor 
merging or star formation) to the mass build-up. 

4.4.2. Central and effective densities 

In Fig. 9 we plot the evolution of the effective and cen- 
tral density in the central galaxies as a function of time. 
The effective density is defined as the mass within the ef- 
fective radius divided by the spherical volume within the 
effective radius. This density is affected by both mass 
and size changes and decreases by 1-2 orders of magni- 
tude between z = 1.5 and z = 0. The density within the 
inner 2 physical kpc, on the other hand, stays roughly 
constant over the last 9 Gyr of cosmic evolution. Major 
mergers mildly increase this central density while mass 
losses by stellar winds tend to decrease it gradually. This 
constancy is not an artifact of our minimal star formation 
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Fig. 7. — The origin of the stellar mass contained within 20 kpc of the central galaxies. The stellar mass fractions are normalized to 1. 
(Left panel) Stellar particles that are found in the central galaxy at a given specified redshift. (Right panel) Stellar particles that have been 
accreted or formed in the specified redshift interval. In each panel the sub-diagrams refer to Gl (left), G2 (middle) and G3 (right). The 
different colors correspond to a different origin of the stellar mass: formed before z = 4 (green), in situ star formation (red), star formation 
in a merging satellite (magenta), formed in a satellite and accreted in a major merger (dark blue) or minor merger (light blue), formed in 
a satellite and accreted after being stripped from the satellite (cyan). The fraction of accreted stars formed in unresolved substructures is 
negligible. 



TABLE 7 

Linear regression parameters of mass-size evolution. 



SFcorr 



Group 


(Mo yr-i) 


2^max 


a 


Gl 





0.54 


1.44±0.05 


G2 





0.40 


2.06±0.03 


G3 





0.69 


1.03±0.03 


G2-HR 





1.6 


1.03±0.02 


Gl 


6 


0.82 


1.56±0.04 


G2 


6 


0.54 


2.09±0.02 


G3 


6 


0.82 


0.97±0.04 


G2-HR 


6 


1.6 


1.04±0.02 


Gl 


11 


1.00 


1.56±0.04 


G2 


11 


0.58 


2.20±0.03 


G3 


11 


0.82 


0.95±0.06 


G2-HR 


11 


1.6 


0.98±0.03 


Gl 


16 


1.00 


1.48±0.05 


G2 


16 


0.64 


2.24±0.04 


G3 


16 


0.79 


0.88±0.09 


G2-HR 


16 


1.6 


0.95±0.04 



Note. — The columns denote: (1) the name of the simulated 
galaxy group, (2) the amount of central star formation subtracted 
according to the minimal correction scheme (a value of means no 
corrections, while we use SFcorr =11 is the default value), (3) the 
highest redshift at which the effective radius is resolved (i.e. i?eff > 
2xebar)5 (4) the best fit exponent of a robust linear regression to the 
mass-size relation from Zmax to z ~ and its error (one standard 
deviation). 



correction scheme since it remains even in the case of no 
correction. We infer that the central density of massive 
galaxies today should correspond closely to their central 
density at z ~ 1.5, while the effective density is strongly 
evolving. 

Between z ~ — 1.5 the mass in the central 2 kpc 
stays roughly constant, but the galaxies increase their 



total stellar mass by a factor of 3-4. Therefore, they 
need to accrete or form mass outside this central region. 
To explore this issue we plot in Fig. 10 the time evo- 
lution of the surface profile of the stellar mass outside 
2 kpc. The increase in the mass surface density of the 
envelope is not smooth, but undergoes phases of fast and 
slow growth. For example the envelope of G2 is almost 
non-evolving between z = 1.5 and ^ = 1, but increases 
rapidely before and after this period. A similar behaviour 
can also be observed at different times for Gl and G3. 
The growing stellar envelope can roughly be fitted with 
a deVaucouleurs profile, although deviations at both low 
radii, e.g. caused by the central "bulge" or the forma- 
tion of a stellar disk, and at large radii, e.g. due to tidal 
debris from satellites, are visible. 

We want to caution the reader that despite the con- 
stancy of the stellar mass Mcen within a small radius 
Teen ("constant central density") the profile within Vcen 
could evolve. For example, one could imagine that vio- 
lent relaxation processes during the initial collapse stage 
and/or by subsequent merger events establish a mass sur- 
face profile of deVaucouleurs form. In this case one can, 
for each given effective radius reff, adjust the r = sur- 
face mass density S(0) such that the mass within Teen is 
Mcen- At fixed Mcen a larger rgfr implies a smaller S(0), 
a larger E at large radii, and an increase in total mass. 
Plugging in Teen = 2 kpc and Mcen ^ 10^^ M0, the model 
predicts that a change in effective radius from 1 kpc to 
5 kpc goes along with an increase in total mass from 
^ 1.5 X lO^^M© to 2.8 X lO^^M© and an increase in 
surface mass density at r = 10 kpc from 2 x lO^M© kpc~^ 
to 2 X 10^ Mq kpc~^. However, our simulations show that 
the fitted effective radii of the envelopes do not strongly 
change between z ^ 2 and z = 0. We thus conclude that 
at 2; ~ 2 the profiles of the progenitor galaxies should 
differ substantially from a single deVaucouleurs law. 
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FlG. 8. — Masses and half-mass radii of the central galaxies in the simulations Gl — G3 and G2 — HR. The left plot shows the stellar 
mass within 20 kpc (top panel) and the half-mass radii (bottom panel) as function of age of the universe (bottom axis) or redshift (top 
axis). The results for simulations Gl (green dot-dashed), G2 (red solid), G2 — HR (magenta solid) and G3 (blue dashed) are indicated. 
Grey lines are measurements of masses and sizes that are not properly resolved. The point with error bars indicates the typical effect at 
t 10 Gyr of varying the residual SFR between 6, 11 (our fiducial value, see Appendix C) and 16 MQyr~-^. The right plot depicts the 
mass-radius evolution of the different galaxies. The symbols correspond to the specific times in left plot. 
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time [ Gyr ] 

Fig. 9. — Evolution of the stellar densities in the simulations 
Gl — G3 and G2 — HR. (Top panel) the effective density peff = 
M(< reff)/V(< Teff) as function of time. The density evolution is 
determined by both the change of the effective radius and the total 
stellar mass of the central galaxy. (Bottom panel) The central 
density pc = M(< 2kpc)/y(< 2 kpc) remains roughly constant 
indicating that the mass within the central 2 kpc of the central 
galaxy does not change much with time. Varying the residual SFR 
correction by ±5 Mq yT~^ results to first order in an overall shift of 
the plots by ~ 0.3 dex, but does not change the general behaviour. 
Symbols are as in Fig. 8. 



4.4.3. The driver of the size growth 

We now address how the stehar envelope is built and 
thus which processes drive the size growth. To this end 
we measure the accreted stellar mass, the stellar mass 
redistributed through the effective radius and the mass 
formed by in situ star formation as function of time. The 
top row in Fig. 11 shows that mass accretion rates in 
mergers can be extremely high (~ lOOOM0yr~^) for a 
very short period of time (~ 10 Myr). This number is 
consistent with what one would expect for groups with 
velocity dispersions of a few hundred km/s, satellite sizes 
of a few kpc and masses of order of lO^^M©. The middle 
row indicates that (i) the accretion flux is due to re- 
solved, singular merger events and not due to a smooth 
accretion of halo stars, (ii) that the accretion by merging 
is predominantly driving the mass evolution and (iii) in 
periods in which no mergers occur either in situ star for- 
mation {z ~ 0.2 — 0.5 in G2) or a redistribution of stellar 
mass {z - 0.7 - 0.8 in Gl, z - 0.1 - 0.2 in G3) can 
drive the mass deposition in the region i?eff-20 kpc and 
hence the size evolution. The bottom row demonstrates 
that, when accumulated over the history of the central 
galaxy, accretion of stars by merging is the dominant 
contributor of stellar mass outside Rei^{z = 0) ~ 3.5 kpc, 
and hence the dominant mechanism that determines the 
sizes of group central galaxies by z = 0. Averaged over 
the three groups merging contributes^ to 70l^5% of the 
mass outside Rei^{z = 0). In addition, non-central star 
formation {14^.1^%) and a redistribution of preformed 

^ The mass fractions are averaged over the simulations G1-G3. 
Quoted as uncertainties are the differences between the average 
mass fraction and the largest and smallest mass fraction among 
the three groups. 
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Fig. 10. — Evolution of the surface mass density profiles of the central galaxies in the groups Gl, G2, G3 (from left to right). A 
deVaucouleurs profile is fit to the data between Sxe^ar cind up to 20 kpc. The effective radii that are derived from the fit are denoted in 
the top right corner of each panel. The profile within Sxe^ar is affected by numerical resolution and this region is shaded gray to caution 
the reader. A particular time span of significant evolution in the surface mass density is shown for each of the central galaxies (blue and 
red symbols). 



stellar mass are also of significance, especially 

in phases without merging activity. For the latter, secu- 
lar, process we observe a correlation between the reshuf- 
fling of mass and the orbiting of satellites in/through the 
central regions of the respective group, indicating that 
satellites may be directly involved in heating the stellar 
component via tidal shocks or by dynamical friction. In 
addition, we also see that the reshuffling is substantial 
for the central galaxy of group G2 when it forms a large 
stellar disk. In this case satellites may induce stellar bars 
or spiral waves and indirectly lead to a redistribution of 
the angular momentum of the stellar component (Sell- 
wood & Binney 2002; Roskar et al. 2008b, a; Kazantzidis 
et al. 2008; Sanchez-Blazquez et al. 2009), see column 3 
of Fig. 2. On the other hand, the flux of redistributed 
stellar mass into the region between 3.5 and 20 kpc (bot- 
tom row of Fig. 11) is smaller by a factor of about 2 in 
G2 — HR compared to G2. This indicates that numerical 
resolution may, at least partially, affect the reshuffling of 
stellar matter, e.g. by changing the bar strength. 

4.5. Beyond mass: The role of assembly history 

The halo masses and radii of the three studied groups 
Gl — G3 below z = 1 differ by less than a factor of 2. 
At z = 0.7 the virial masses are (by selection) almost 
identical, while at z = they range from 1.1 to 1.6 x 
IO^^Mq, see Fig D4. Nonetheless, we observe strong 
differences in morphology of the central galaxies at z < 1, 
see Fig. 2 and Fig. 3. In this section we will discuss 
the connection between the morphology of the central 
galaxies and their mass assembly histories and cooling 
properties. 

4.5.1. Star formation histories and the cold gas 

reservoir 

In order to understand the differences between the 
three central group galaxies we begin by investigating 
their star formation histories. The star formation his- 
tories peak at z ~ 4, when the universe was about 2 
Gyr old (left panel of Fig 12). Later several short bursts 
of star formation occur of which some are connected to 
merging events. The star formation history of group G2 
differs significantly from the groups Gl and G3. The 
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Fig. 11. — Mass transfer into a shell outside i^eff- Different rows 
correspond to the different groups: Gl (left column), G2 (middle 
column, thin lines), G2 — HR (middle column, thick lines), G3 
(right column). (Top row) the net fluxes in the region between R^^ 
and 20 kpc averaged over 13 Myr; (Middle row) as the top row but 
averaged over 100 Myr; (Bottom row) cumulative net fluxes into 
the region between Reff(z = 0) ~ 3.5 kpc and 20 kpc. The different 
lines correspond to the net fluxes of accreted stars (blue solid line), 
of stars formed in situ outside the effective radius (green dash- 
dotted line) and of stars formed in situ within the effective radius 
but redistributed outwards (red dashed line), respectively. 



overall star formation rate in the former group is gener- 
ally larger, especially at high redshift (z ~ 3 — 4) and be- 
low z = 1. At high redshifts the star formation timescale 
{M^/SFR) of the three groups at a given instant of time 
is similar. It ranges from ~ 0.9 Gyr at z = 3 (i.e. small 
compared to the age of the Universe of 2.2 Gyr at z = 3) 
to 2 — 5 Gyr at z = 2 (which is similar or even longer than 
the age of the Universe of 3.4 Gyr at z = 2) - indicat- 
ing that the importance of star formation reduces with 
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time and becomes subdominant around z = 2. Despite 
their similar star formation timescales the central galax- 
ies reach a stellar mass of IO^^Mq at different times: 
G2 at 2: ~ 2.4, while Gl and G3 cross this boundary 
much later at around z ~ 1.5. When we compare the 
star formation history with the available cold gas mass 
within 20 kpc we generally see a co-evolution (compare 
left and middle panel of Fig 12). The groups cross the 
predicted virial mass threshold of Mghock ~ 7 x 10^^ for 
the generation of a stable shock near the virial radius 
(Birnboim & Dekel 2003) at z - 2.3 (Gl), z - 2.9 (G2) 
and z ~ 2.0 (03). Consistent with the transition from 
a cold to a hot gas accretion mode (Dekel & Birnboim 
2006, Keres et al. 2005, 2009, Ocvirk et al. 2008) we find 
that the cold gas mass within the central 20 kpc of the 
groups drops roughly at these predicted times, as seen 
in the middle panel of Fig 12. The time delays of 0.5-1 
Gyr between the decrease in available cold gas and the 
reduction in the star formation rates compare well with 
the gas consumption timescales Mcoidgas/SFR ~ 0.2 — 0.4 
Gyr and indicate that the star formation rates mainly de- 
crease due to the lack of cold gas in the central galaxies. 
Between z ^ 1 and z ^ 0.2 group G2 contains about 
1 X IO^^Mq of cold gas within 20 kpc. This exceeds the 
cold gas reservoir of Gl and G3 within that time-span 
by an order of magnitude and explains why the central 
galaxy in G2 experiences a large star formation activ- 
ity at those intermediate redshifts of about 40-50 Mq 
yr~^. By z ~ 0.5 it resembles a galaxy with prominent 
gaseous and stellar disks with spiral features and with 
a red bulge-like component. A bulge-to-disk decompo- 
sition with GALFIT results in B/D-^ 1-1.7, depend- 
ing on projection and whether mass density or I-band 
images are fitted. At the same time the rotational sup- 
port of the galaxy smoothly increases from v/cr = 0.3 at 
z ^ 1.1 to V / a = 1.1 dit z ^ 0.5. On the other hand the 
central galaxies in Gl and G3 are relatively quiet (SFR 
< IOM0 yr~^) during that periode and their star forma- 
tion is restricted to the central ~ 1 kpc. Below z ^ 0.2 
the cold gas in G2 depletes quickly and its star formation 
rate drops substantially. Simultaneously, minor mergers 
and tidal interactions with smaller satellites in the group 
heat and damage the stellar disk and convert the spiral 
galaxy into an SO over a timescale of 1-2 Gyr. Hence, the 
simulations indicate that, firstly, the cold gas fractions 
drop and thus the formation of new stars is suppressed, 
before subsequently major galaxy mergers (Gl, G3) or 
minor mergers and tidal interactions with satellites (G2) 
heat the stellar system and transform it towards a more 
early-type galaxy. 

4.5.2. The origin of the cold gas 

The substantial amount of cold gas of G2 below z = 1 
originates from three different sources: (i) cooling from 
the hot gas halo, (ii) cold gas brought in by mergers 
or in small gas clumps (but not by cold streams), and 
(iii) cooling from a warm, primordial gas phase. The 
first process contributes the most. More than 70% of 
the cold gas at z = 0.5 has cooled out of the hot halo 
below z = 1. The cooling time of the hot gas within 
the central 10 (20) kpc of group G2 is - 100 (300) Myr 
below z = 1, see Fig. D6. Therefore cooling can pro- 
ceed quickly unless heating mechanisms keep the gas hot, 
such as adiabatic compression due to a smooth accre- 



tion of intergalactic gas onto the halo, shock heating of 
the ambient gas due to supersonic collisions of infalling 
gas lumps (e.g. Ryu et al. 2003), heating by dynami- 
cal friction of supersonic galaxy motions (Ostriker 1999; 
El-Zant et al. 2004), clumpy accretion (Dekel & Birn- 
boim 2008), heating by a background UV field (Haardt 
& Madau 1996; Hoeft et al. 2006), or heating due to 
feedback mechanisms such as radio mode AGN (Croton 
et al. 2006). The latter is not included in our simula- 
tions. Roughly 15% of the cold gas reservoir cools from 
a warm (~ 10^ K), zero-metallicity (primordial) gas that 
is accreted from the surroundings of the group halo, pos- 
sibly from the left-overs of cold filaments (Keres & Hern- 
quist 2009) or from underdense, unprocessed gas. The 
accretion of metal enriched, cold (T < 3.2 x 10^ K) gas 
contributes to ~ 10%. At z ^ 0.4 group G2 undergoes 
a merger with another, smaller group of virial mass of 
~ 5 X 10^^ which has not yet crossed the threshold mass 
above which its cold accretion is terminated. It carries a 
cold gas mass of a few 10^^ Mq. However, most of the 
cold gas in that subgroup is shock heated or consumed by 
star formation during the group-group merger, so that it 
barely contributes to the cold gas reservoir of the main 
halo. 

4.5.3. The cooling time and assembly history 

We now address why the cooling in the group G2 and 
in the groups Gl and G3 behaves so differently. We 
believe that the mass assembly history, in particular the 
frequency of major halo-halo mergers, impacts the gas 
cooling in at least two different ways: 

First, at fixed halo mass the absence of late halo-halo 
mergers implies an earlier formation time and therefore, 
on average, a more concentrated mass profile of the halo 
(Navarro et al. 1997; Wechsler et al. 2002). The more 
concentrated halo has a higher central density and may 
be more easily prone to cooling instabilities. In dark- 
matter-only re-simulations of the groups G2 and G3 the 
concentration^ of group G2 amounts to about twice the 
concentration of group G3, namely c = 5.9 vs c = 2.8, 
see Fig. D5. It is straightforward to estimate that for 
an NFW profile a doubling in concentration increases 
the enclosed mass at small radii by a factor 2-2.8 (for 
c ^ 2 — 20). However, this is clearly an oversimplifi- 
cation given the non-linear properties of the gas cooling, 
the adiabatic contraction of the dark matter and the fact 
that baryons and dark matter might decouple in the cen- 
tral region. In the right panel of Fig. D6 we plot den- 
sity weighted temperature, density, entropy and cooling 
time of the gas within 20 kpc for our three groups. At 
2; ~ 0.7 the central gas density of group G2 excels the 
central densities in the groups Gl and G3 by a factor 
~ 4, the central temperature of G2 is smaller by 10-20%, 
and consequently its cooling time is lower (300 Myr vs. 
> 1 Gyr). 

The second impact of halo-halo mergers it that they 
contribute to the heating of the intra-group gas. For 
example, infalling satellites and gas clumps can convert 
their potential and kinetic energies into thermal energy of 

^ The concentration c is defined as the ratio between virial radius 
(for the calculation of c we use the definition of Wechsler et al. 
(2002)) and scale radius of the best fitting NFW-profile (Navarro 
et al. 1997). 
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time[Gyr] time[Gyr] time[Gyr] 

Fig. 12. — (Left panel) The star formation rate within 20 kpc around the central galaxy of the groups Gl (green dot-dashed line), G2 
(red solid line), G2 — HR (magenta solid line) and (blue dashed line) a function of cosmic time and redshift. The star formation is 
corrected for an artifical central star formation rate of ~ IIMq yr~^, see appendix C. (Middle panel) The mass of cold gas within the 
inner 20 kpc of the central group galaxies. (Right panel) The dark matter mass within a fixed physical radius of 200 kpc around the main 
progenitor. 



the ambient medium (Johansson et al. 2009). This heat- 
ing mechanism is boosted in a halo-halo mergers that 
bring in a large number of satellites. But also other 
channels such as large scale shock heating of the gas 
may be important or the decrease of the gas density by 
outward motions of the gas (increase in kinetic energy). 
In the right panel of Fig. 12 we plot the dark matter 
mass within 200 kpc around the central galaxy. Com- 
pared with the virial mass as function of time (Fig. D4) 
this plot allows to identify both the first pericenter of 
the merging halos and the final coalescence, leading to a 
characteristic "U" shape. Together with Fig. D6 (right 
panel) this plot shows that the halo-halo mergers that 
occur between z = lA and 2: = 1 in both Gl and G3 
lead to an increase in temperature and decrease in den- 
sity at the centers of the groups. Similarly, the merger 
that occurs between z ^ 0.5 and z = 0.25 in group G3 
terminates a short cooling event at z ~ 0.5, and the 
merger that starts at z ~ 0.3 in group G2 terminates the 
cooling below z ~ 0.25. 

The mass accretion history is therefore an important 
factor in order to explain the diversity that group halos 
of the same mass show in their cooling properties or the 
related properties (cold gas fractions, star formation ac- 
tivity, morphology, etc.) of their central group galaxies. 
On the one hand, halo-halo mergers inject energy into 
the intra-group gas and prevent or even terminate hot 
accretion episodes. On the other hand, groups with late 
halo-halo mergers, with a higher substructure fraction or 
with later formation times have on average less concen- 
trated dark matter halos (Navarro et al. 1997; Wechsler 
et al. 2002; Gao et al. 2004), which increases their cooling 
times. We thus anticipate for groups of a given mass on 
average an anti-correlation between the amount of sub- 
structure (related to formation time) and the star forma- 
tion rate of the central galaxy. On the same grounds we 
expect that groups which undergo mergers less likely host 
a disk or SO galaxy at their centers compared to well viri- 
alized groups with a quiet merger history. We note that 
the assembly bias for the formation time and concentra- 
tion is small at the group mass scale (Gao & White 2007; 
Li et al. 2008). Therefore, the presented mechanism does 



not strongly differentiate between groups in underdense 
regions and groups near a cluster, respectively. We note 
that our analysis does not include the impact of AGN 
feedback on the galaxy group scale, which is an impor- 
tant, yet open question for future work. 

5. SUMMARY 

We have simulated the evolution of galaxy groups with 
a final virial mass of ~ IO^^Mq in order to study the 
evolution of their central galaxies. We ran the Tree-SPH 
code GASOLINE (Wadsley et al. 2004) with the same 
parameters and at a resolution comparable to the one 
used in current state-of-the-art simulations that follow 
the formation and evolution of disk galaxies in a cos- 
mological context (Governato et al. 2007). All our cen- 
tral group galaxies end up with roughly the same stellar 
mass of a few times 10^^ Mq at z = and their stel- 
lar profile can be well fitted with a Sersic index of > 4. 
Furthermore, our simulated central galaxies are kinemat- 
ically hot systems supported by velocity dispersion, have 
colors of evolved stellar populations and an early-type 
morphology. However, the detailed morphology (E or 
SO) and the amount of rotation differs among the cen- 
tral galaxies. We trace such differences back to their 
different merger and gas accretion histories. Overall the 
basic properties of our simulated central galaxies match 
approximatively that of observed central group galaxies 
in the local Universe (Yang et al. 2008; Guo et al. 2009; 
Cibinel et al. 2009a,b,c). 

At redshift ~ 2 the most massive progenitor galaxies 
reside in dark matter halos of galactic size (~ lO^^M©). 
The parent halos are just massive enough to support a 
stable shock at the virial radius that reduces the im- 
pact of smooth cold accretion (Dekel & Birnboim 2006). 
However, some cold streams are still able to reach the 
central objects and so at z ~ 2 the galaxies experience 
still a significant (> 20 Mq yr~^) star formation that is 
sustained by a substantial reservoir of cold gas replen- 
ished by cold and hot accretion modes. At this time the 
galaxies have gathered, within a compact (~ 1 kpc), blue 
(unless it is heavily dust obscured) "proto-bulge" , an al- 
ready significant fraction (10-25%) of their final stellar 
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mass. They can be observed as K^^^^ ~ 20 star forming 
BzK galaxies or (as long as Ay < 1) as UV-optical se- 
lected galaxies (Daddi et al. 2004; Adelberger et al. 2004; 
Grazian et al. 2007). Only in case they are strongly dust 
obscured {Ay ^ 1.3) they would qualify at z > 2.3 as 
distant red galaxies (Franx et al. 2003). The cold gas 
arranges in form of extended gas disk of a few kpc ra- 
dius in agreement to observations of massive, star form- 
ing disks (Genzel et al. 2006, 2008; Kriek et al. 2008a; 
Forster Schreiber et al. 2009). Over the next two Gyr 
the galaxies gradually reduce their amount of star for- 
mation and the overall colors redden due to the aging 
stellar population. The galaxies are still very compact - 
they concentrate ~ 10^^ Mq within 1 kpc at z ~ 1.5 and 
might thus contribute to compact galaxies observed at 
those redshifts. In the remaining 9 Gyr of cosmic time 
the galaxies are subject to major and minor mergers, 
episodes of gas accretion and a secular redistribution of 
pre-formed stellar mass that drive the growth of a stellar 
envelope around a now rather passively evolving "bulge" . 

The merger histories of the three central galaxies and 
their progenitors differ substantially. One galaxy grows 
its stellar mass since z ^ 1.5 mainly by minor merging 
and star formation. In the other two groups the cen- 
tral galaxies experience each two major mergers below 
redshift z ^ 1. Two out of the 4 occurring major merg- 
ers involve dissipative disk-disk mergers at z ~ 0.8 and 
z ^ 1 and produce gas-poor, pressure supported rem- 
nants which subsequently undergo a "dry merger" with 
another gas-poor, pressure supported galaxy at z = 0.4 
and z = 0.1, respectively. 

Despite the diversity in their merger histories the evo- 
lution of the stellar component of the progenitors of the 
central group galaxies proceeds in a similar fashion. The 
densities within the inner 2 kpc remain roughly constant 
since z ^ 1.5 and the galaxies subsequently grow in size 
by building a stellar envelope around the early formed 
stellar bulge. The growth of this envelope reduces the 
density within the half-mass radius by almost 2 orders 
of magnitude since z ^ 1.5. It is interesting to com- 
pare our result with the proposed size evolution of com- 
pact, passively evolving, red-and-dead galaxies observed 
at z ~ 1 — 2 (e.g. Cimatti et al. 2008; van Dokkum et al. 
2008). In particular, minor merging has been advocated 
as a driver since it is expected to lead to a stronger size 
evolution per unit mass than major merging (Naab et al. 
2009; Bezanson et al. 2009) and is potentially required 
to connect the sizes of compact, passive galaxies at high 
z to the local early-type galaxy population (Bezanson 
et al. 2009). Although minor mergers happen numer- 
ously in our simulations the size growth is overall slower 
than predicted from minor merging alone. The simple 
reason is that major mergers occur frequently (in our 
case in 2 our of 3 simulated galaxy groups) and then 
naturally dominate the mass evolution. Also, in one of 
our groups non-central star formation temporarily drives 
a late mass and size growth. We expect this latter chan- 
nel to become dominant for lower mass galaxies (see e.g. 
Guo & White 2008). In addition, we expect central star 
formation to slow down the size growth. Our results im- 
ply that the size growth per unit mass of massive, yet 
star forming galaxies at z ~ 2 is qualitatively similar, 
but likely less steep, compared to the mass-size evolu- 
tion of massive, compact and passively evolving galaxies 



observed at ^ ~ 2. 

Below z = 1 the main progenitors of the central group 
galaxies differ substantially in their morphologies and 
star formation activities although they sit in groups of 
similar virial mass (cf. Scannapieco et al. (2009) for 
galactic halos). In one of three cases the central galaxy 
accretes a substantial amount of gas between z ~ 0.3—0.9 
that cools out of the hot intra-group medium due to the 
short cooling time at the center of the group. This leads 
to an inside-out growth of a stellar disk and the mor- 
phology approaches that of a massive spiral galaxy. The 
cooling episode is eventually stopped when a halo merger 
with an infalling subgroup occurs. The central galaxy is 
subsequently transformed into an early-type (SO) galaxy 
as cooling from the hot gas halo is suppressed, the cold 
gas reservoir is exhausted by star formation, and mi- 
nor mergers and tidal interactions with satellites dam- 
age and heat the stellar disk below z = 0.2. In con- 
trast, the other two studied groups are subject to several 
group-group mergers which provide an energy injection 
mechanism that increases the temperature and decreases 
the density of the intra-group gas and thus counteracts 
a hot accretion mode. In addition, group mergers are 
followed by major galaxy mergers between the former 
central galaxies of the subgroups creating remnants of el- 
liptical morphology. Hence, in the absence of additional 
heating mechanism the mass accretion history plays an 
important role in limiting the strength of the hot accre- 
tion mode and thus the kinematics and morphology of 
the central group galaxies. If the cooling is suppressed 
galaxy mergers and tidal interactions with satellites are 
able to drive a rejuvenated central galaxy towards an 
early-type appearance. 

6. OPEN ISSUES 

A missing ingredient in our simulations is feedback 
from active galactic nuclei typically associated with 
super-massive black holes which are found in the centers 
of most galaxies (Kormendy & Richstone 1995; Magor- 
rian et al. 1998; Gebhardt et al. 2000; Ferrarese & Merritt 
2000). AGN feedback could potentially rival with envi- 
ronmental effects and assembly history in setting the final 
galaxy properties, although recent work shows that this 
is not necessarily the case (Croft et al. 2009). Indeed our 
results argue against a major role of AGN feedback, at 
least for the properties of central galaxies outside the in- 
ner kpc. This is because even without AGN feedback our 
simulations reproduce reasonably well many properties 
of z ^ early-type galaxies (cf. Khalatyan et al. 2008), 
such as the colors, morphologies, kinematics and struc- 
ture outside the ~ central softening lengthy albeit with 
masses and/or sizes that are somewhat biased w.r.t. av- 
erage values seen in local surveys (Shen et al. 2003). The 
tendency towards compactness of our simulated galaxies 
is the counterpart of the mass concentration problem still 
seen in the simulations of disk galaxies, whereby, despite 
the fact that disk scale lengths and overall disk sizes are 
reasonable reproduced at the highest resolution achieved 
so far, central bulges are still too massive and compact, 
see Mayer et al. (2008). This likely reflects the need for 
even higher resolution in the early stages of galaxy as- 
sembly, which are responsible for setting the central den- 
sity and inner mass distribution, as shown clearly in this 
paper, and/or the necessity of additional heating mech- 
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anisms in the very innermost regions. Therefore, the 
study of AGN feedback and its effect on massive galaxies 
in galaxy groups may be the key to obtain more realis- 
tic central densities, and thus effective radii, by partially 
suppressing cooling and star formation in the central kilo- 
parsec. Alternatively, the same effect may be achieved 
by a more realistic modeling of the multi-phase interstel- 
lar medium and of star formation/supernovae feedback 
once molecular gas densities are resolved (Robertson & 
Kravtsov 2008, Governato et al. 2009b). 

Our results suggest that halo-halo mergers are an effi- 
cient means to quench cooling. The precise mechanism 
(or the chain of mechanisms) by which the potential en- 
ergy of the pre-merged halos is converted into kinetic or 
thermal energy of the gas still needs to be worked out, 
and we do not exclude that major stellar mergers are 
involved. However, we want to point out that a corre- 
lation but not a causal connection between shutdown of 
hot accretion and major galaxy merging may arise even if 
major mergers are not responsible for the suppression of 
the cooling - simply because after the halo-halo merger 
terminates the cooling, the former central galaxies will 
engage in a major merger due to the short dynamical 
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friction time for massive galaxies. 

Another question that remains open is how the (al- 
ready massive) galaxies at 2; ~ 2 relate to higher redshift 
objects, such as 2; ~ 3 Lyman-Break galaxies (Steidel 
et al. 1996) or Lyman-Alpha emitters (e.g. Hu et al. 
1998) that are a potential progenitor population given 
their clustering properties (e.g. Blanc et al. 2008). We 
intend to address this question with future simulations 
that are specifically designed to resolve the z > 3 pro- 
genitors of central group galaxies. 
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APPENDIX 

STAR FORMATION CORRECTION SCHEMES 

The raw star formation rates of the central galaxies are about ~ 11 Mq yr~^ at z = and are therefore an order of 
magnitude higher than the average star formation rates of massive galaxies in the local Universe (Zheng et al. 2007). 
This high star formation biases masses, sizes, colors and luminosities of the studied galaxies. The star formation excess 
can (i) arise from a shortcoming of the implemented model, e.g. a lack of a physical heating mechanism such as radio 
mode AGN feedback (Croton et al. 2006), (ii) be caused by a too coarse resolution to resolve the relevant heating 
mechanism, such as clumpy accretion (Dekel & Birnboim 2008) or (iii) be of purely numerical origin. Unfortunately 
we do not know the origin of the excessive star formation in our simulations, but we can study its properties and 
try to correct for it in a posteriori fashion. Ultimately, however, the simulation model and/or resolution should be 
improved in order to reduce the high star formation rate in a fully self-consistent manner. In Table Al we summarize 
the consequences of various star formation correction schemes on masses and sizes of the central galaxies. 

As Fig. A2 demonstrates a significant amount of star formation is occuring within the unresolved centers (r < ebar) 
of the galaxies. While within a radius of 2 x ebar the amount of SF varies substantially with time, group and resolution 
the SFR within r < Cbar depends only little on resolution and studied group. It still varies with time although it 
hardly drops below 11 Mq yr~^. A natural way to correct for excessive central star formation is therefore to either (i) 
subtract a constant star formation rate of the order of 11 ± 5 Mq yr~^, or (ii) to subtract the complete star formation 
within ebar- In both approaches we can correct stellar masses by subtracting an excess- mass Mexc(^) from the central 
region that arises from central star formation since z = 4 and takes mass losses by stellar winds into account. 

Mexc(^) = / dt' SFR{t') uo{t - t') 

Here, uj{t) denotes the mass of a star particle at time t after its formation (see Fig. Al). If we assume a constant 
star formation rate of 11 ± 5 Mq yr~^ for 2: < 4 we obtain an excess mass of 8.4 ± 3.8 x IO^^Mq. Half-mass radii are 
computed after the excess mass is substracted from the inner region (r < 2 x ebar)- 

We call the first approach the minimal star formation correction approach since its effects are relatively weak: the 
masses, luminosities and sizes at z = change by less than 25%. The colors (Table 4) are hardly affected by any star 
formation correction method because, as described in section 3, colors are effectively measured at a radius of > 2 x ebar 
by fixing the mass-to-light ratio of the central, unresolved region to its value at 2 x ebar- 

The second approach has a stronger impact on SFR, masses (and sizes), and luminosities. Comparing Fig. A3 and 
Fig. 12, for example, demonstrates that the SFR is strongly affected by the second correction scheme. In particular, 
central galaxies are essentially evolving completely passive below z = 2 except for short bursts associated with mergers 
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Fig. Al. — The mass fraction of a star particle, representing a simple stellar population with Scalo IMF, that is removed due to stellar 
winds vs. its age. The mass loss reaches 40% after several Gyr. 

and a few continues, long lasting gas accretion episodes. The in situ formed stellar mass reduces by a factor of ~ 3 
compared to the no-correction scheme and by a factor of ~ 2 compared to the minimal correction scheme, see Fig. 
A4. The correction is somewhat more important for the in situ formed stellar mass of Gl and G3 and less for G2. 
However, the central galaxies in Gl and G3 build their mass mostly by merging. We do not know how the amount of 
artificial star formation in less massive galaxies/halos and we can only speculate on how a change in the in situ formed 
mass is refiected on the total stellar mass. If we simply correct for the amount of in situ formed stellar mass within 
ebar of the central galaxy, the stellar mass will reduce by only 30 — 35%. If we assume, however, that the in situ star 
formation in all progenitors, not only the main progenitor, needs the same correction factor, the masses are reduced 
strongly by a factor ~ 3. A complication of the second correction scheme is that it is not clear how galaxy sizes have 
to be corrected. 

We chose to employ the minimal correction approach partially for this reason, partially in order to avoid overcor- 
rections and also because the results we focus on (e.g. the exponent in the mass-size relation) do not seem to strongly 
depend on which correction we apply. However, the absolute masses and sizes of the produced z = galaxies are 
in much better agreement with observations if a strong star formation correction is applied (c.f. the discussion in 
section 4.2, Table Al). We caution the reader that our correction scheme neglects the dynamical effects of having an 
excess mass in the center and that of artificially enhanced gas infiow from larger radii. The excess mass changes the 
potential of the galaxy, decreases the dynamical time of orbiting matter, and the cooling time of gas in the hot diffuse 
medium near the center. The excess mass is thus likely to increase non-linearly over time. In addition, artificially 
enhanced central gas cooling can lead to a stronger baryonic contraction of the stellar and dark matter and hence to a 
reduction of the half-mass radii of the central galaxies. Better resolution (especially at higher redshifts), a higher star 
formation threshold or an even more efficient feedback in the lower mass progenitor galaxies could obliterate the need 
for an explicit correction of the star formation and help to remedy any remaining discrepancies between simulated and 
observed galaxy properties. 

METAL COOLING 

High-temperature metal cooling is expected to modify the star formation rates and masses of galaxies by increasing 
the rate at which gas can cool. Our simulations of galaxy groups are actually at a rather sweet spot (Tvir ~ 10'' K) in 
which the cooling rate is only mildly dependent on metallicity (Sutherland & Dopita 1993). In contrast, the cooling 
rate differs by more than an order of magnitude at the virial temperatures typical of galactic halos. Quantifying the 
impact of metal cooling is a non-trivial issue, because higher star formation leads to enhanced stellar feedback that 
will inject more energy into the ambient gas. Numerical simulations that study this non-linear interplay between gas 
cooling, accretion and feedback find that the inclusion of high-temperature metal cooling increases the cosmic star 
formation rate by 0.25-0.3 dex at 2; ~ 1 and 0.1-0.16 dex at z ~ 3 (Choi & Nagamine 2009; Schaye et al. 2009), and is 
hardly significant at higher redshifts (Hernquist & Springel 2003; Schaye et al. 2009). While metal cooling increases 
the stellar masses of low and intermediate mass galaxies, it leaves the stellar masses of the most massive (> 10^^ Mq) 
galaxies at 2; ~ 1 almost unchanged (Choi & Nagamine 2009). This is partially due to the lower metal-sensitivity 
of the cooling function at high temperatures and also because more massive galaxies form on average their stars at 
earlier times when the impact of metal cooling is reduced. Metal cooling extends the cooling function to temperatures 
far below 10^ K. We do not expect that this affects the masses of the galaxies, but it might have an effect on central 
concentrations as gas is allowed to get colder and denser. However, the Jeans mass of the 300 K cold gas phase is small: 
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Fig. A2. — The star formation rate within 1 (sohd hnes) and 2 (dotted lines) baryonic softening lengths of the simulations Gl (green), G2 
(red), G2 — HR (magenta) and G3 (blue). The horizontal dashed line indicates the estimate ~ 11 ± 5Mq yT~^ of residual star formation 
rate in the central region. 
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Fig. A3. — The star formation rates of the central galaxies in the groups Gl (green), G2 (red), G2 — HR (magenta) and G3 (blue) within 
a spherical volume of 20 kpc around the central region after subtracting the star formation occurring within a baryonic softening length. 



~ 5 X 10^ Mq at an unrealistically low density of 0.1 cm~^ and 10^ — 10^ Mq at the typical densities (~ 200 cm~^) 
of a giant molecular cloud. Even our high resolution simulation is far from resolving these Jeans masses and thus not 
able to study reliably the properties of such a cold gas phase. We note that, for consistency, our star formation and 
stellar feedback parameters are constrained by disk galaxy simulations with primordial gas composition (Stinson et al. 
2006) and chosen such as to reproduce the observed Kennicutt-Schmidt law. 

RESOLUTION TESTS 

Throughout the paper we present the results of simulation G2 alongside with the results obtained from its higher 
resolution counterpart G2 — HR. In this section we present additional considerations regarding the sizes and masses 
of the simulated central galaxies. 

The measurement of masses and sizes is potentially affected by, e.g. overcooling, heating due to two-body relaxation, 
by the softening of the gravitational potential or by artificially enhanced mass deposition in the central region. In Fig. 
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Fig. A4. — The in situ formed stellar mass within a radius of 20 kpc around the center of the central galaxy. (Top panel) raw simulation 
output, (Middle panel) after correcting for an artifical star formation rate of 11 Mq yr"-*^, (Bottom panel) after correcting for artificial star 
formation within the inner baryonic softening length. Mass losses by stellar winds are taken into account in the mass correction. 



TABLE Al 

Masses and Sizes of the central group galaxies at z=0 with and without star formation corrections. 









^eff 


Group 


SFcorr 


(IQiiMo) 


(kpc) 


Gl 


no 


4.8 


3.1 




subtract 11 ± SM© yr~i 


4.0 ±0.4 


4.0 ±0.5 




subtract < Cbar (central) 


3.2 






subtract < Cbar (all) 


1.6 




G2 


no 


5.3 


2.7 




subtract 11 ± SM© yr~i 


4.5 ±0.4 


3.7 ±0.5 




subtract < Cbar (central) 


3.4 






subtract < Cbar (a^ll) 


1.8 




G3 


no 


5.1 


2.5 




subtract 11 ± SM© yr~^ 


4.3 ±0.3 


3.2 ±0.4 




subtract < Cbar (central) 


3.7 






subtract < Cbar {six) 


1.7 





Note. — The first column denotes the studied galaxy group. The second column states the correction method for artificial star formation 
that is applied: The top row (no) refers to no correction. The next row refers to the minimal correction scheme. The third row shows the 
result of a scheme that removes in situ star formation of the central galaxy (and its most massive progenitor) within one baryonic softening 
length. The 4th row refers to a correction scheme in which the total stellar mass of the central galaxy is corrected with the same factor 
as its in situ formed stellar matter, see text. Columns 3 and 4 denote the stellar mass within 20 physical kpc and the three-dimensional 
radius that contains half of Mtot- 

CI we therefore compare the stehar masses produced at intermediate (G2) and at high resolution {G2 — HR). Below 
z = 2 the stellar mass M^, the stellar mass within rmin ~ 2 x ebar and the stellar mass outside rmin (and also the 
spherical and projected mass density profiles outside 3 x ebar, not shown) converge to within a ~ 10 — 20% accuracy. 
Hence, we require that the half-mass radius is (at least) of the size ~ 2 x ebar for a reliable size measurement. In order 
to avoid that size determinations depend too strongly on the mass within rmin we also employ another safety criterion. 
We estimate the mass within rmin by fixing the density at rmin and extrapolate the density inwards with a power law of 
slope -1. This density profile corresponds to a fiat surface density and we find that under these corrections the masses 
within rmin and half-mass radii change both by less than 30% at z = 1. However, at z = 1.5, the mass is affected at 
the 50% level and the half- mass radius by a factor of 2. We therefore restrict the study of half-mass radii to z < 1.5. 
In Table CI we compare the structural and kinematic properties of the central galaxies in G2 and G2 — HR at z = 0.1 
(the lowest redshift reached in the simulation G2 — HR) finding good agreement. 
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Fig. CI. — Comparison of the stellar masses of the central galaxies between the intermediate resolution G2 and its high resolution 
counterpart G2 — HR. The black solid line shows the mass difference of G2 and G2 — HR normalized to the mass of G2 — HR (all within 
20 kpc). The red dot-dashed line shows the fractional difference of stellar mass within the inner 1.44 kpc, i.e. the unresolved region of G2, 
while the blue dashed curve shows the fractional mass difference between 1.44 kpc and 20 kpc. 



TABLE CI 

Structural and kinematic properties of the central group galaxies at z ■ 



: 0.1. 



Group (IO^'^Mq) 



Rvir 

(kpc) 



(kpc) 



odeVau 
^eff 



oSersic 
^eff 



(kpc) (kpc) n 



(km/s) 



Vrot/cTefi 



G2 11.8 420 41.9 ±3.7 3.1 ±0.5 3.9±0.3 4.4±0.3 3 332±20 0.87±0.05 

G2-HR 11.4 416 39.6 ± 3.6 2.9 ± 0.4 3.2±0.2 2.7±0.2 5 348±20 0.92±0.05 

Note. — The columns are as in Table 3. The kinematical and structural properties of the central galaxy in the simulations G2 and 
G2 — HR are both inferred from radii that exclude the central 2-3 baryonic softening lengths in units of the simulation G2. 




Fig. D1. — The number density of galaxies exceeding a given stellar mass derived from the stellar mass function of Baldry et al. (2004). 
In addition to the integrated number density of all galaxies in the local universe (solid black line) the contributions from red (dashed red 
line) and blue (dot-dashed blue line) galaxies according to the color-magnitude split of Baldry et al. (2004) are shown. The exponential 
cut-off at the high mass end causes the number density to drop by an order of magnitude when the stellar mass doubles from 2 x lO-*^-*^ Mq 
{n = 1.4 X 10-4) to 4 X 10^1 Mq {n = 1.8 x lO"^). 
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Fig. D2. — Apparent Js — Ks color vs. Kg magnitude of the main progenitors of the central group galaxies within a projected radius 
of 8 kpc between z = ?> and z = 1.5. The different lines corresponds to Gl (green), G2 (red), G3 (blue) and G2 — HR (magenta). The 
horizontal, dashed lines indicates the color cut that separates distant red galaxies (DRGs; Jg — Kg > 2.3) from non-DRG (Js — Kg < 2.3). 
Gray symbols indicate the change of the colors due to different extinction corrections (Calzetti et al. 2000) with Ay ranging from to 2. 
The color evolution is shown for the default value Ay = 0.8. Colors and magnitudes are normalized to the Vega system using the stellar 
template of Kurucz (1992). 
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Fig. D3. — Apparent G — IZ vs. Un — G colors of the main progenitors of the central group galaxies within a projected radius of 8 kpc 
between z = 3 and z = 1.5. The different lines corresponds to Gl (green), G2 (red), G3 (blue) and G2 — HR (magenta). The selection 
windows of the BM and BX techniques are indicated as colors boxes. Gray symbols indicate the change of the colors due to different 
extinction corrections (Calzetti et al. 2000) with Ay ranging from to 2. The color evolution is shown for the default value Ay = 0.8. 
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Fig. D4. — The virial mass of the groups Gl (green dot-dashed), G2 (red sohd), G3 (blue dashed) and G2 — HR (magenta sohd) as 
function of time (bottom axis) and redshift (top axis). 
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Fig. D5. — The density profiles of the SPH groups Gl (green dot-dashed), G2 (red solid), G?> (blue dashed) and G2 — HR (magenta 
solid) and the density profiles of corresponding dark-matter-only re-simulations: G2 — DM (black solid), G3 — DM (cyan dashed). The 
profiles of G2 - DM and G?> - DM are best fitted by NFW profiles (colored dotted lines) with c = 5.9 (G2) and c = 2.8 (GS). The fit 
ranges are indicated by vertical dotted lines. 
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Fig. D6. — (Left) Radial profile of density-weighted temperature, number density of charged particles, entropy and the cooling time of 
the intra-group medium in the group G2 — HR at different times. Temperature, density, entropy and cooling time are calculated from 
either all gas particles (thin dashed line) or from only the hot gas (> 3.2 x IC^ K) particles (thick solid line) in spherical shells around the 
group center. (From top to bottom): (1) The density- averaged temperature initially raises quickly and by z ~ 2 has reached an almost 
constant value of ~ 10^ in the central region, while it continues to raise in the outer region as the halo grows. The hot accretion below 
z ~ 1 is clearly visible as a drop in temperature within the inner 10 kpc, (2) The number density of charged particles (electrons, protons 
and helium ions) n = p/ji under the assumption of complete ionization of a zero-metallicity gas {ji = 0.59 amu), (3) The entropy of the 
gas S = Tj'n?!'^, (4) The cooling time (eq. 14 of Scannapieco &; Oh 2004) increases with increasing distance and becomes longer than the 
age of the Universe (dotted horizontal lines) at a radius of 100-200 kpc. (Right) Temperature, density, entropy and cooling time (defined 
as in the left panel) calculated from either all gas particles (light gray lines) or from only the hot gas (> 3.2 x 10^ K) particles (colored 
lines) within a sphere of 20 kpc around the centers of the groups G\ (green dot-dashed), G2 (red solid), 6*3 (blue dashed) and G2 — HR 
(magenta solid) as function of time (bottom axis) and redshift (top axis). In all three groups the cooling time within 20 kpc is always lower 
than the age of the Universe (dotted line). 



